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Abstract. We present numerical simulations of radiative 
transfer in the spiral galaxy NGC 6946. The interstel- 
lar medium is represented as a two phases medium, with 
molecular clouds and a smooth diffuse phase. The molec- 
ular gas distribution is calculated in a self-consistent way 
from the distribution of an ensemble of molecular clouds 
evolving in the gravitational potential of NGC 6946. We 
simulate star formation by creating OB associations in 
molecular clouds. The transfer of UV radiation is calcu- 
lated in the clumpy interstellar medium, to determine the 
local UV illumination of molecular clouds. We compute 
the emergent intensity in the UV continuum (912 - 2000 
A), in the Ha and 0+ ^^3/2 — ^ -P1/2 lines as well as in the 
continuum at far infrared wavelengths, 60, 100 & 200 /xm. 

It is possible to obtain a consistent picture of this 
galaxy with a global star formation rate of 4 MQyr"^ (for 
stars with masses in the range 2-60 Mq) occuring mostly 
in the spiral arms. The close spatial association of massive 
stars and molecular clouds has a profound impact on the 
transfer of UV radiation in the galactic disk and on the 
dust emission. The median distance travelled by UV pho- 
tons is about 120 pc. However, when they have escaped 
from the vicinity of their parent OB associations, UV pho- 
tons may travel quite far in the disk, up to 1 kpc. The UV 
opacity of the model spiral galaxy disk, observed face-on, 
is 0.8 at 1000 A and 0.7 at 2000 A. 

For radii less than 4 kpc, the C^ 158 fim line is mostly 
produced in photodissociation regions at the surfaces of 
molecular clouds. The C"'' emission from diffuse atomic 
gas accounts for about 20% of the total. It becomes sig- 
nificant at large distance from the nucleus (r > 4 kpc). 
Molecular clouds and diffuse atomic gas have almost equal 
contributions to the total far infrared emission from 60 to 
200 fim. As a whole, 72% of the 60 - 200 FIR emis- 
sion can be attributed to dust grains heated by the UV 
radiation of massive stars and 28% by the radiation field 



of the old stellar population. 
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1. Introduction 

Because the peak of the thermal dust emission occurs 
in the far-infrared (60-200 fim) for Galactic interstellar 
clouds, this wavelength domain is recognized as a very im- 
portant 'window' to study the interstellar medium, both 
in the Milky Way and in external galaxies. However, due 
to the lack of spatial resolution provided by telescopes 
operating at these wavelengths, the interpretation of the 
signals measured in external galaxies is still a subject of 
controversy. The debate is focussed on the origins of the 
heating of dust grains and of the ^^3/2 — ^ P1/2 C+ line at 
158 fim. This fine structure line of ionized carbon is one of 
the main cooling lines of neutral atomic gas (Wolfire et al. 



1995), but it is also seen in HII regions, with different exci- 
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tation mechanisms in these two cases. Other sources of C+ 
emission are the dense photodissociation regions (PDRs) 
where intense UV radiation from massive stars impinges 
on molecular clouds. In these regions, the photodissocia- 
tion of carbon monoxide creates a layer of ionized carbon 
at the edge of the molecular cloud, which is a source of 
intense C^ emission. 

The far infrared emission is produced by dust grains 
heated by star light but the relative contribution of young 
massive stars on one hand, and of the bulk of the stel- 
lar population on the other h and, a re still debated (see 
e.g. W alterbos & Gr eenaw alt 1996 , Dever eux fc Young 
1990|, Co x & Mezger |l989i Thronson et al. |1990| , Rice et 
al. 1990 ). Using low resolution data from IRAS, IR col- 
ors have been used to assess the respective roles of star 
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forming regions and cirrus in external galaxies at a global 



scale (e.g. Calzetti et al. 1995). It is difficult to find spi- 
ral galaxies dominated by a single population of heating 
sources, since the relative contributions of massive stars 
and of the disk population vary from galaxy to galaxy, 
and probably from place to place in a given galaxy. A 
main problem of these studies is the lack of spatial resolu- 
tion at far infrared wavelengths. Even when it is possible 
to use high spatial resolution data at millimeter and sub- 
millimeter wavelengths, the gap in spatial resolution from 
the local interstellar clouds to external galaxies is huge: 
for nearby galaxies located at distances D of a few Mpc, 
a one arcmin beam encompasses 0.3(D/lMpc)kpc, much 
larger than the size of a single molecular complex in the 
Galaxy, whereas structures down to 0.01 pc are commonly 
observed in local clouds (Falgarone et al. 1992) and in the 
Magellanic Clouds (Rubio et al. |1993| ). 

Numerical simulations provide a way out of this prob- 
lem: it is possible to reproduce the observed molecular 
gas distribution of nearby spiral galaxies with numerical 
simulations using observed data as input parameters. For 
example Garcia-Burillo et al. ( 1993| ) were able to fit the 
spatial distribution and kinematics of the molecular gas in 
M 51 usi ng the cloud collision code developped by Combes 
& Gerin ( 1985 ). Though this code has a spatial resolution 
of a few hundred parsecs (the cell size for the large scale 
dynamics), it is possible to include "micro-physics" at the 
parsec scale inside each cell. We have taken this approach 
and implemented star formation in this code to study the 
far infrared and C"'" emissions in the spiral galaxy NGC 
6946. The next section summarizes the current data on 
NGC 6946. We describe the model in Sect. 3 and present 
the results for NGC 6946 in Sect. 4. The implications of 
this work for the interpretation of the C+ and FIR emis- 
sion of spiral galaxies are discussed in Sect. 5. 

2. NGC 6946 

NGC 6946 is a nearby Scd galaxy with a low inclina- 
tion, i = 34° (Considere & Athanassoula 198^ ). At the 
adopted distance of 5 Mpc, 1 arcmin on the sky corre- 
sponds to a linear size of 1.5 kpc. The main characteris- 
tics of NGC 6946 are summarized in Table 1. It shows a 
weak bar (Martin 1995) and a prominent spiral structure 
with both an m = 2 and an m = 4 pattern (Considere & 
Athanassoula 19881 ). With this open spiral structure, the 
arm and interarm regions are resolved at moderate spatial 
resolution (30" ) . The gas distribution is well known since 
NGC 6946 has been mapped at high spatial resolution in 
HI (B oulanger & Via llefon d |l99|), and CO (Casoh et al. 



giola 1991). It also belongs to the normal galaxy sample 
studied by ISO and as such has been extensively mapped 



1990| , Clausset et al. |l99lD . Figure 1 presents a C0(2-l) 



map obtained at the IRAM 30m telescope (Clausset et al. 
1991), and Fig. 2 a V-I color map obtained at the Ob- 



servatoire de Haute-Provence (OHP) (P. Boisse, private 
communication). NGC 6946 is marginally resolved in the 
far infrared maps obtained by IRAS and the KAO (Engar- 



in the mid and far infrared (Malhotra et al. 1996, Helou 
et al. p^96|, Lu et al. p^96|. Tuffs et al. |T996|). NGC 6946 is 



forming stars actively, as revealed by its bright Ha emis- 



sion (Bonnarel et al. 1988a, 1988b, Kenn icutt 1989) and 



t he nu merous supernovae (Li & Li 1995 ). Madden et al. 
(1993) have mapped the C+ 158 fj,m emission at 45" reso- 
lution with the KAO. Bright 0+ emission is found in the 
nucleus and in the disk. At this scale, the C+ emission is 
correlated with the CO emission in the spiral arms and 
with the HI emission in the outer disk. 



3. A model of NGC 6946 

3.1. Representation of the interstellar medium. 

We have chosen to include two phases for the interstellar 
medium : 

i) a dense molecular phase composed of spherical molec- 
ular clouds with masses ranging from 10^ to 10^ M© and 
mean density 50 H2 cm~^. This mean density is used to 
determine cloud sizes from their masses. The half power 
width of the molecular layer, defined as < >^/^, is equal 
to 65 pc. The global volume filling factor of this molec- 
ular phase in the disk, (H2), amounts to 2.5% . The 
clumpy structure of molecular clouds is taken into account 
by using a larger density, 5 x 10^ H2 cm^^, to compute the 
C+ emission. The 0+ emission is not sensitive to the gas 
density when it is l arger than ^ 10'^ cm~^ (see Tielens 



Hollenbach 1985 and Sect. 3.6). Indeed, C+ observa- 



tions of nearby molecular clouds show that C+ emission 
is detected over t he who le projected surface of molecular 
clouds (Jaffe et al. 1994 ). The total mass of dense molecu- 
lar gas in the model is set to match the value deduced from 
CO(l-O) observations with a conversion factor of 2.3 x 10^° 
H2Cm-2 (K kms-i)-i, 2.0 x 10^ Mq (Table 1). 
ii) a neutral diffuse phase, with a nearly constant HI den- 
sity. We have adjusted the value of the mean HI density us- 
ing the data from Boulanger & Viallefond ( |1992| ): it varies 
slowly from 0.8 H cm""^ at the center to 2 H cm~'^ in the 
outer disk {r 6 kpc). We have assumed that the HI disk 
has a vertical thickness of 360 pc — 2Hz, and that the 
density is uniform from —H^ to +Hz. The velocity dis- 
persion of the diffuse gas is set to 10 kms^^. The total 
mass of diffuse neutral gas inside a radius of 6 arcmin is 
taken from the HI observations to be 2.6 x 10^ M©. 

In addition to these two phases, we take into account 
the formation of HII regions around OB associations (see 
below). For the numerical calculations, we represent the 
galaxy disk as a large 3D grid. The total grid length L 
and the cell size rceii are chosen to match the following 
constraints: 

- the grid size has to be at least as large as the optical 
disk, which has a diameter of 12.5 kpc, in order to be able 
to follow the path of photons in the external regions of the 
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galaxy, and the heating of the neutral gas at large distance 
from the nucleus. 

- the spatial resolution, rceiu should be good enough to 
describe the clumpy structure of the medium: in partic- 
ular the smallest clouds in the ensemble should occupy 
at least one cell. With the adopted parameters for the 
clouds, the radius of the smallest clouds with a mass of 
10^ Mq is 4.5 pc. This number is only an estimate of the 
actual size of a given cloud. It is possible to use an alter- 
native method based on the scaling relations for molecular 
clouds, the so-called Larson's laws. Using the m ass-ra dius 
relation M(M0) = IQQBj^ {^c^) (Falgaroneet al. |l992D , the 
radius of a 10'^ M© cloud is then 3 pc. Another constraint 
on the spatial resolution is provided by the size of OB as- 
sociations, which is generally larger than 30 pc (Garmany 



1994) and their distance from molecular clouds. Leisawitz 



(1991) has estimated the mean distance between OB as- 
sociations and their parent cloud to be about 50 parsecs. 

We have chosen as the best compromise to represent 
the galaxy with a 1024 x 1024 x 24 cell structure. The 
galaxy size is then L = 2 x Rmax = 12.5 kpc and the 
resolution Vceii — 12.2 pc. The cells are filled with neu- 
tral atomic gas, or with molecular gas at the molecular 
cloud positions, an 10^ Mq molecular cloud then occu- 
pies 3x3x3 cells. Finally, we include the ionized gas 
in the Stromgren spheres centered on each OB associa- 
tion, and replace the molecular and atomic gas with ion- 
ized gas whenever required. The radii of the Stromgren 
spheres are calculated assuming classical ionization in HI- 
bounded HII regions (Miller & Cox 1993). Because of the 
coarse spatial resolution of the dynamical code (200 pc), 
we do not attempt to reproduce the interstellar medium 
in the central region of NGC 6946 (r < 500 pc). 

3. 2. Spatial distribution of molecular clouds 

This section of the code uses the cloud-cloud collision code 



described by Combes & Gerin (1985), and used by Garcia- 
Burillo et al. ( |1993D and Gerin et al. ( |l99lD to model the 
gas dynamics in nearby galaxies. The molecular clouds 
move in the gravitational potential of the galaxy, they 
grow through cloud-cloud collisions, which are sticky pro- 
cesses, and are disrupted by simulated SN events. The 
molecular gas is immediately recycled into small molec- 
ular clouds. The gravitational potential is deduced from 
an R band image of NGC 6946 (Viallefond & Bonnarel, 
private communication, Bonnarel et al. 198§| ), with fore- 
ground stars removed. The image has been rotated to put 
the major axis vertical, and deprojected to face-on, using 
as projection parameters PA = 69° and i = 34 °. The 
image has then been binned to 256 x 256 pixels. This R 
band image covers 8.5' x 8.5' on the sky, corresponding 
to a radius of 6.24 kpc at the assumed distance of NGC 
6946, 5 Mpc. The actual spatial resolution amounts to a 
few times the cell size, about 200 pc. The following step 
is to build a good gravitational potential from this image. 



As in Garcia-Burillo et al. ( 1993| ), the calculation is done 
in two steps. The axisymmetric part of the potential is 
obtained by assuming a constant mass to light ratio, and 
adjusting this constant to reproduce the observed CO and 
HI rotation curves. The perturbations due to the bar and 
spiral arms are included in the non axisymmetric part of 
the potential. The last parameter to adjust is VLp, the pat- 
tern speed of the density wave (bar -I- spiral arms) . The gas 
distribution and velocity field are very sensitive to Sip. We 
find that Q,p= 42 kms^^ kpc^^ gives the best results. The 
corotation is found at a radius of 3.5 kpc, close to the end 
of the bar, and the OLR lies outside of the disk. Figure 3a 
presents the deprojected R-band image. Figures 3b and 3c 
present an example of the obtained molecular cloud distri- 
bution. The gas clouds are concentrated in the spiral arms, 
with few molecular clouds at radii larger than 4 kpc. The 
model rotation curve is shown in Fig. 4, together with the 
angular frequency ^{r). We have not attempted to model 
the compression of diffuse gas in the spiral structure, and 
have kept an axisymmetric distribution of diffuse atomic 
gas. 

3.3. Stellar population 

Once the gas is distributed over the galactic disk, we cre- 
ate OB associations, with star masses ranging from 10 to 
60 Mq (stellar types from B2 to 05). We have chosen not 
to include stars less massive than 10 Mq because their life- 
time becomes a significant fraction of the rotation period. 
Shortward of 2000 A, the UV radiation is mostly produced 
by OB stars and the contribution by late B type stars is 
at most 20% (Waltcrbos & Greenawalt 1996, Mathis et al. 



1983). The choice of the upper mass cutoff at 60 Mq is mo- 



tivated by the work of Heydari-Malayeri & Beuzit (1994) 
who have shown that suspected very massive stars (~ 100 
Mq) are actually clusters of less massive stars. Also, a 
very short time step is required to sample adequately the 
lifetime of these very massive stars. 

We allow only clouds more massive than 4 x 10^ Mq 
to form massive stars. In the cells satisfying this criterion, 
the OB associations are born at the outer edge of the 
molecular clouds. We have tested different star formation 
laws, so as to match as well as possible the Ha radial pro- 
files. The best fits are obtained with a star formation rate 
(SFR) depending on the angular frequency and the local 
gas surface density as proposed by Wyse & Silk (198£) : 



SFR{r)^en{r)agas{r) 

with e the star formation efficiency, ^{r) the angular fre- 
quency at radius r and Ugasij") the HI+H2 surface density 
at position r in the disk. We find that e — 5% provides the 
best match to the Ha radial profile by Kennicutt ( 1989| ). 

The stellar mass distribution inside an OB association 
is drawn from the Initial Mass Function (IMF). We have 
chosen an index close to the Salpeter IMF, as suggested 
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by measurements in Galactic and extragalactic OB asso- 
ciations (Massey et al. |l99|) : ^^^^ = Af-^-^. 

The stellar associations are born with a mean veloc- 
ity relative to their parent cloud chosen from a Gaussian 
of mean value 10 kms-i and FWHM ~ 10 kms'^ These 
figures are in good agreement with observed data by Lei- 



sawitz et al. (1989) in their study of the relation of star 
clusters with molecular clouds. With this value, the aver- 
age motion of OB stars relative to their parent cloud is 10 
pc in 10^ years. Finally, stars die after a time equal t o thei r 
Main Sequence lifetimes from Glisten & Mezger ( 1983 ). 
Effective temperatures, stellar luminosities and radii, av- 
eraged over the main-sequence lifetime, are taken from 
Cox et al. ( |1986 ). Because we are only interested in broad 
band fluxes and colors, we have calculated all stellar fluxes, 
from UV to B bands, in the black-body approximation. 
Lyman continuum radiation production rates, averaged 
on the Main Sequence lifetime, are also taken from Cox 



et al. (1986) and Giisten & Mezger (1983) to have a co- 
herent set of parameters. We use the parametrization as 
a function of the stellar mass, M: 



;o5io(^) = 38.3 



.16 X logw{M) 



-0.24 X logl^^{M) - 0.41 x logl„{M). 

to compute the thermal radio continuum flux. Ha lumi- 
nosity and Stromgren sphere diameters. 

We also include the effects of the radiation of massive 
stars on neutral gas, and allow molecular clouds to be par- 
tially eroded and ionized by the radiation of nearby OB 
associations, if they overlap with the Stromgren sphere of 
an OB association. We compute the radius of this sphere, 
Rhii , assuming that all stars are located at the same po- 
sition and have a global production rate of Lyman contin- 
uum photons N^yc, equal to the sum of the contribution 
of the individual stars, and for case B recombination in a 



diffuse medium of density uh- Rhii = ( 



3 NLy 



4 TTQ 



)i/3. In 



this formula, as is the hydrogen recombination coefficient 
(Osterbrock |1989| ). The mean radius for Khii is 35 pc. 

To constrain the population of massive stars, we calcu- 
late different stellar outputs, namely we perform a detailed 
calculation of the UV radiation field at 912-2000 A, and 
also compute the global U flux and U radial profile, as 
well as the thermal radio continuum and Ha emission. 

• 6 cm radio-continuum. 



According to Mezger ( |1972| ) and Turner & Ho ( |1994| ), the 
radio continuum emission of HII regions for case B recom- 
bination, at an electronic temperature Tg of 10'* K and 
with 45 % of the ionizing photons being converted into 
Ha photons, is directly related to N^yc by: 

= 1.1 X 1050(^)(- ^ ^2 



' mjy Mpc 

This formula does not include any correction for dust ab- 
sorption of the ionizing radiation within the HII region. 



Current estimations are that nearly 50% of the Lyman 
continuum may be absorbed by dust. 
• Ha. 

The relationship between the Ha luminosity and the pro- 
duction rate of Lyman continuum photons, using the same 



assumptions as above, can be deduced from Mezger (1972) 
and Peimbert et al. (1975): 



L 



Ho 



= 3.55x 10-46(^). 



The correction for the extinction may amount to about 1 
magnitud e at Ha but it is highly uncer tain (Mc Kee & 
WiUiams 1997, van der Hulst et al. 1988). Because of the 



uncertainties involved in this correction, we have preferred 
not to do it. This allows us to check that the energy is 
conserved with a good accuracy in the simulation. 

The formed OB associations have typical production 
rates of Lyman continuum photons, N^yc, in the range 
, The cumulative Lyman continuum lumi- 
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nosity function of the population is shown in Fig. 5. Its 
shape is similar to th e dist ribution for Galactic HII regions 
(McKee & Williams 1997). There is a small excess around 



^Lyc > 10^^-^ s~*, dominated by the most massive stars 
in the younger associations, and a deficiency of associa- 
tions with N^yc > 10^°-^ s~*. As a whole, the match of 
the two distributions is very good. The distribution of the 
intrinsic UV luminosity of OB associations extends over 
four orders of magnitude from lO'' to lO'^ Lq, i.e. from 
small associations with about 20 stars and a total stellar 
mass of 600 Mq, up to large associations gathering 200 
OB stars, including a few 50-60 Mq stars, and having a 
total stellar mass of 6000 M©. 



3.4- Dust properties 

We use the averag e Gal actic extinction curve from Fitz- 
patrick & Massa ( 1988 ) at UV wavelengths. For visible 
and near infrared wavelengths, we use the work by Seaton 
( 1979 ). This curve is probably valid on a large scale in the 
diffuse medium of NGC 6946 which has a similar metal- 
licity to the Milky Way. We assume a constant metallicity 
and gas to dust ratio in the disk, and use the average value 
for the Milky Way : e(s^ = 5.8 x lO^i H cm'^ 
and Ay — 3.1 Eb-v (Bohlin ct al 



1978) 



mag 
The extinction 
through a cell filled with molecular gas, with the adopted 
spatial resolution, is 2.5 magnitudes. The dust properties 
have been summarized by Bruzual et al. ( 1988| ) and Witt 
& Gordon ( 1996| ). Apart from the enhanced absorption in 
the 2175 A bump, the dust albedo is fairly constant at UV 
and visible wavelengths at a; ~ 0.55. We include the coher- 
ent scattering of UV light by dust grains. The anisotropy 
is described using the Henyey-Greenstein (1941) function, 
where the anisotropy parameter is defined as g —< cos9 > 
and 6 is the scattering angle. 
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3.5. Radiative transfer 

The observed properties of external galaxies depend on 
the propagation of the stellar radiation in the interstellar 
medium. To determine the local radiation field, we fol- 
low the propagation of UV photons (912-2000 A) emit- 
ted from the OB stars in the two phases medium. In 
each cell, we compute a local radiation field. We define 
the local radiation field in the UV, xuvi relative to the 
mean radiation field in the UV at the Solar radius, xoj 
estabhshed by Mathis et al. ( |1983| ), so that xuv = 



In the Galaxy and at the solar radius, the InterStellar 
Radiation Field, ISRF in the UV, has been defined as 



Xo 



4ttJ = 4tt f^™?-^ JxdX = 1.84 x 10 



J9I2A 



erg cm ^s 

/•2/im 



•'912A 



and through visible/IR bands as Gq 

2 X 10~^ erg cm~^s~^, where is the specific intensity 
of the radiation field averaged over Att sr (Mathis et al. 



1983) 



On their path, photons can be scattered and/or ab- 
sorbed by dust in both neutral phases. We have chosen 
not to follow individual photons which would have been 
time consuming, but to use instead pseudo-photons rep- 
resenting a collection of N photons. We are then able to 
probe a larger area with a lower number of photons. We 
sample the 912-2000 A interval with 20 bins of constant 
wavelength width AA, and launch n^ pseudo-photons per 
association per wavelength bin. The longward limit has 
been set to 2000 A, to avoid having to take into acount 
the contribution from star types later than B to the ISRF 
(Walterbos & Greenawalt, 1996). 

The pseudo-photons emitted by the OB association 
number i carry a fraction fl dX of the luminosity L\ radi- 
ated by this OB association. They travel from the center 
of a cell (size rceii) to an adjacent one, and the energy 
absorbed by the interstellar medium, when non-zero, is 
left on the common face of these 2 cells. The absorbed en- 
ergy is reprocessed in the far-infrared. The pseudo-photon 
energy along its path from association number i can be 
written as : 

fldX^a'^ 

In this formula, a* = (x, y, z) is the product of 

the probabilities for non absorption in each cell along the 
travel from the OB association number i to the cell po- 
sition (x,y,z), or with uj(X) the dust albedo and tx the 
total cell opacity including scattering : 
Q,i = g-(i-'^(>'))TA 

in the diffuse phase, and 

= a; (A) 
in the H2 phase. 

When leaving an OB association, the direction for each 
pseudo-photon is uniformly chosen on the unit sphere. 
These pseudo-photons travel in the two-phase interstellar 
medium. There are different possibilities when reaching a 
new cell : 



i) The gas in the cell is diffuse and atomic. Then the 
pseudo-photon can either : 

- be absorbed partially in the cell. A fraction of the lu- 
minosity is left, and the pseudo-photon continues in the 
same direction with a lower luminosity. 

- be scattered. 

ii) The cell is filled with molecular gas. 

The pseudo-photon is partly absorbed, and partly back- 
ward scattered, the ratio between the energy left in the cell 
and the total energy of the pseudo-photon depends on the 
albedo as l — uj{X). We have chosen backward scattering in 
that case to prevent the pseudo-photon from interacting 
with the same molecular cloud several times. 
After a scattering event, a scattering angle 9 is chosen 
according to the anisotropy function, and the aziniuthal 
angle 4> is uniformly chosen in the interval [— tt, tt]. Then 
the new direction is easily deduced from the previous one 
(Witt |1977D . 

Since we cannot store all the pseudo-photons' incident 
directions, and in order to define an isotropic radiation 
field in each cell, we assume that the equivalent surface 
of the cell is 6r^g;;. Then, the UV radiation field in the 
cell (x, y, z) due to the contributions of all pseudo-photons 
travelling through this cell, can be expressed relative to 
the Galactic ISRF xo as : 



xix,y,z) 
Xo 



i SX 



Because of the coarse angular resolution and the low num- 
ber of pseudo-photons leaving each OB association, the 
resultant map of the radiation field presents strong fluc- 
tuations. We have chosen to smooth the map of the ra- 
diation field by averaging the data in nearby cells, typi- 
cally 3x3 cells. Furthermore, some cells are never visited 
by UV photons, for example in the inter arm region or in 
the outer disk (r > 5 kpc). In that case, we use as incident 
radiation field, the Galactic ISRF (longward of 2000 A), 
scaled by the local surface brightness in the R band image 
to take into account the radial variation of the radiation 
field from the old stellar population between the central 
regions (r < 1 kpc) and the outer disk (r > 5 kpc). The 
reference value is given in units of 47r T^^™. JxdX — Gn - 

Xo — 1-82 X 10^^ erg cm^^s^^. The position for the ref- 
erence value has been chosen at the edge of a spiral arm, 
at a distance r = 4 kpc from the center, where the UV 
radiation field xuv is close to 1. 

3.6. Emergent emission. 

Once the local UV energy density has been calculated, 
models are used to determine the C+ and FIR emissions 
from each cell, assuming that they arise from the same 
area as the one used for the calculation of the UV radiation 
fiel d. For the dust emission, we use the model by Desert et 
al. ( 1990| ), which has 3 different components: PAHs which 
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are fully ionized when xuv > 1, very small grains and big 
grains. We calculate the dust emission in the four IRAS 
bands at 12, 25, 60 and 100 plus an additional band 
at 200 fj,m, as the reprocessing of the combination of the 
UV radiation field (912 - 2000 A) described by xuv and 
of the ISRF for the 2000 A-2 /xm part of the spectrum. 
We decrease the PAH abundance in large r adiat ion field 
environments, as suggested by Ryter et al. ( 1987 ), to one 
fifth of the standard value when xuv > 100. 

Very large molecular clouds, with masses larger than 
10^ M0 occupy more than 4 cells in the grid. The inner 
cell is not directly exposed to the UV radiation and for 
this cell we assume that the ISRF is attenuated by 2.5 
magnitudes of visual extinction. 

We assume low optical depth in the mid and far 
infrared. From the analysis of the COBE maps of the 
Galaxy, Boulangeret al. (|1996D deduce tx/Nh= 1 x lO'^s 



cm (A /250 /J,m)~ , which combined with the mean 

column density of individual clouds, = 1.5 x 10 
cm^^, gives an opacity of 2 x 10"'^ at 200 fim. 

According to the dominant phase in a given cell of the 
model, the emergent infrared emission can be : 

- reprocessing of the whole incident stellar radiation for 
molecular clouds which are totally optically thick in the 
UV. All impinging radiation is completely reprocessed in 
infrared emission from the outer cells of molecular clouds, 

- proportional to the gas column density for the diffuse 
medium, which is optically thin in the UV. We assume no 
UV extinction at the 12 pc scale. 

We use the PhotoDissociation Region (PDR) model by 

Pl/2 C + 



Le Bourlot et al. (1993) to estimate the R 



3/2 ■ 



emission at the surfaces of molecular clouds. We use a con- 
stant molecular hydrogen density of 5 x 10'^ H2 cm~"^ and 
the incident UV field xuv ■ The density is not a critical pa- 
rameter as long as it is higher than the critical density for 
collisional excitation of C"*" (1000 H cm~'^) (see Tielens & 
HoUenbach 1985). Furthermore, we assume that the whole 



surface of clouds contributes to the 0+ emission. As for 
the C"*" emission from the diffuse neutral phase, we use 



the model by Wolfire et al. (1995) for a two phase neutral 



atomic medium. Finally, the contribution from the ionized 
gas in the HII regions around the OB associations is also 
included. We assume that the gas has the same density as 
the diffuse medium, rie = Uh and an electronic tempera- 
ture of 10"* K. The total C"*" luminosity from an HII region 
of radius Rq is then proportional to the volume of the HII 
region with a correction factor to take into account the 
other ionization stages of carbon. 

Because we deal with a line, the opacity may not be 
small depending on the gas distribution and viewing ge- 
ometry. In fact, opacity effects are important at large in- 
clination angles. To obtain an edge-on view of the model 
galaxy in the ^P3/2 — ^ P1/2 C"*" line, we have made an ac- 
curate calculation of the radiative transfer in this line. For 
each line of sight through the disk, we sample the line pro- 
file with bins of 1 kms^^ width, and calculate the emergent 



intensity in each velocity bin, including saturation effects. 
We assume that the intrinsic velocity dispersion of a PDR 
is 1 kms~^. For the diffuse medium, the velocity dispersion 
is 10 kms~^. We do not account for absorption in the far 
infrared continuum. We obtain tc+ = 0.40 for an edge-on 
view, and Tp+= 0.10 for a face-on view. 

3. 7. Implementation of the model 

To avoid transient stages of the simulation, the code is 
evolved during a few time steps. We stop the simulation 
when stable results are obtained on a time scale of 20 
Myrs. This time scale corresponds to about half the life- 
time of a giant molecular cloud before disruption by pho- 
toevaporation. This is the reason why we can not integrate 
further in time without treating gas recycling. The time 
step has been fixed at 10^ years, shorter than the lifetime 
of the most massive stars. We have checked the reliability 
of the calculations by different tests: 

- We have verified that the total UV luminosity from the 
stellar population emerges from the galaxy either at the 
same wavelength, or at far infrared wavelengths for the 
light reprocessed by dust grains. The total luminosity is 
conserved with an accuracy of 1 %. 

- When the number of pseudo-photons leaving each OB 
association n^ is too small, the map of UV radiation field 
is noisy with a few extremely bright spots and large voids. 
This is due to undersampling of the galaxy volume. The 
number of pseudo-photons should be as large as possible, 
but we have verified that we obtain a good map of the 
UV radiation field with 100 pseudo-photons per OB as- 
sociation. The map is smooth in the vicinity of the OB 
associations, hence the ratio of FIR emissions from the 
diffuse and dense gas stays constant with increasing n^. 

- The cell size is also a critical parameter: since xuv 
is proportional to it might be underestimated for 
small clouds very close to OB associations. This has severe 
consequences for the C"*" emission, which scales roughly 

''^ceii^'^sixuv), but little or none for the FIR emission 
which varies as r'^eiiXuv since in that case there is no 
resultant scaling with Vceii- To test the validity of the 
adopted resolution, we have performed a run restricted to 
one quadrant only, with a cell size of 6.1 pc. We observed 
no large variation in the emission and thus conclude 
that the adopted resolution of 12 pc is correct for our pur- 
pose. Note that the volume filling factor decreases to 1.3 
% in the high resolution run, because we fill the space in 
a more accurate way using a higher spatial resolution. 

4. Results 

Table 2 summarizes the input parameters for the Stan- 
dard Model, and Table 3 presents the results. The star 
formation rate from 2 to 60 M0 is fixed at 4.0 M0yr~-'^, 
with a star formation efficiency e of 5%, as defined in Sect. 
3.4. With these values, the modelled Ha luminosity and 
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Ha radial profile are in quite good agreement with the 
observed data (Kennicutt |1989 ). This is also true for the 
UV luminosity at 2000 A and the 6 cm luminosity. We 
are thus confident that the massive star population is well 
constrained by the observed data. With a low mass cut- 
off at 10 M0 stars, we overestimate the UV luminosity 
at 2000 A of the modelled stellar population by 20%, be- 
cause we miss the contribution to the UV continuum of 
lower mass stars, between 2 and 10 Mq. 



4.I. Disk opacity in the UV 

We have computed an average opacity over the galaxy in 
the UV and for the Ha line. We define this opacity as : 

T = _^„(2,e™er3e«t /^emitted) ^^^^^ j^ermtted ^^le total 

luminosity in the disk at a given wavelength and L^'^'^'^ae-nt 
is the emergent luminosity. This opacity is computed for 
two different viewing angles of the model, z = 0° for face- 
on and i = 90° for edge-on. We have found a significant 
opacity for the face-on view, at 1000 A, 2000 A and for Ha, 
namely t(IOOOA) = 0.8, r(2000l) = 0.7, T{Ha) = 0.60 
for the whole galaxy. 

The opacity is controlled simultaneously by the geom- 
etry of the molecular cloud ensemble and by the diffuse 
medium. If we ignore the extinction due to the diffuse com- 
ponent, we find an opacity of 0.51 at 1000 A. This value 
is due to geometrical effects, mostly blocking of the UV 
radiation by molecular clouds, and it does not depend on 
wavelength. Thus we can write the opacity at any wave- 
length in the UV as t\ = 0.51 + rj^^ , the second term 
accounting for the wavelength dependence of the extinc- 
tion in the diffuse medium. 

A global opacity of t ~ 0.8 corresponds to a fraction of 
approximatively 45% of the far UV stellar radiation leav- 
ing the galaxy disk, mostly above or below the main plane. 
Most of these photons have not been scattered because the 
probability of leaving the disc after a scattering event is 
low. This significant fraction of the radiation from mas- 
sive stars leaking out of HII regions could contribute to 
the maintenance of the Reynolds layer of ionized gas. The 
derived face-on opacity at 2000 A, 0.7, falls well within 



the range of opacities derived by Buat & Xu (1996). The 
mean extinction in their sample of nearby spiral galaxies 
is ~ 0.9 mag at 2000 A. Though the opacity is not very 
large, the mean distance travelled by a UV photon before 
absorption is quite small, 440 pc, roughly equal to the HI 
disc thickness. As shown on Fig. 6, there are however pho- 
tons travelling to much larger distances, 1 to 2 kpc, with 
a small probability (0.01). Conversely, many zones in the 
interarm receive very few UV photons. Due to the lower 
gas density, few OB associations are created in the inter- 
arm region. The numerous OB associations in the arms 
are too distant to contribute to the local radiation field 
since the arm/interarm separation is larger than 1 kpc in 
the disk. 



The distribution of xuv values provide further infor- 
mation on the radiation field resulting from the OB asso- 
ciations (Fig. 7). Whereas most of the galaxy is exposed 
to a low UV radiation field, it is possible to find regions 
with high UV intensity {xuv > 1000) even at a moderate 
spatial resolution. The total dynamical range of the UV 
radiation field extends over more than 4 orders of magni- 
tude. This huge variation can be explained by the close 
association of OB associations and molecular clouds: in a 
galaxy with a prominent spiral structure, OB associations 
are born in the spiral arms, where the gas density is the 
highest. This maximizes both the illumination of molec- 
ular clouds by UV radiation and the absorption of UV 
radiation by molecular gas, hence the heating of molecu- 
lar gas. For the model galaxy, we find that 30 % of the 
total number of cells with molecular gas are exposed to a 
strong or median radiation field [xuv > 10)- These cells 
are located in 40% of the molecular clouds. This figure is 
comparable to the clouds in Milky Way: Solomon et al. 
(1985) found that in the Galaxy, at a resolution greater 
than 10 pc, 25 % of the molecular clouds are warm and 
associated with HII regions. Also, Williams and Mc Kee 



(1997) estimate that at least one OB star is present in half 
of the giant molecular clouds with masses larger than 10^ 
M0 . The probability to find massive stars or clusters asso- 
ciated with a giant molecular cloud increases sharply with 
the cloud mass and reaches almost 1 for m asses larger than 
8 X 10^ M© (Wilhams & MacKee [l997|) . Our numerical 
results are in agreement with these facts. 



4-.2. Far infrared emission 

We now discuss the emergent radiation from the model 
galaxy and start with the infrared emission. As in the 
Desert et al. (1990) dust model, the luminosities in the 
IRAS bands are computed as iirD'^i'Siy, where Si. is the 
total observed flux density and D is the distance to the 
object. The infared colors are given as the flux density 
ratios, to compare with observed data. 

The model galaxy has very similar emissions as 
NGC 6946 at 60-100 & 200 ^im, with outputs of 5.1, 8.5 
and 4.9 x 10^ Lg, corresponding to 114%, 128% and 144% 
of the luminosities observed at those wavelengths. The 
far infrared emission comes from both the molecular and 
atomic gas phases. 

The UV radiation is the main heating mechanism of 
the dense and diffuse gas phases, with contributions of 4.2 
X 10^ L© and 6.1 x 10^ Lq at 60 & 100 /im. The contri- 
bution to the FIR emission of the old stellar population, 
described by the ISRF, is a factor 3 lower, with 0.9 x 10^ 
Lq and 2.4 x 10^ Lq in the 60 & 100 iim bands. The sit- 
uation is different at 200 iim, where dust grains heated by 
the UV radiation or by the ISRF have almost equal contri- 
butions to the total luminosity: 3.0 x 10® Lq for the UV 
and 1.9 X 10® Lq for the ISRF. The contribution from the 
inner parts of clouds illuminated by the attenuated ISRF 
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is only 0.5 x 10* L©. As a whole, 72 % of the far infrared 
luminosity can be attributed to UV heated gas, which is 
mostly molecular. The remaining 28% corresponds to dust 
heated by the ISRF, at locations far away from the OB 
associations. 

The diffuse and dense phases have similar contribu- 
tions to the total FIR emission, with a slight excess from 
the molecular clouds, 54 % versus 46 % from the diffuse 
gas. This significant contribution from the diffuse gas is 
due to the fact that it occupies a large fraction of the 



galaxy volume. Hippelein et al. (1996) also conclude from 
ISO observations of other nearby galaxies (M51, MlOl) 
that the neutral atomic gas has an important contribu- 
tion to the far infrared emission. The contribution from 
the atomic gas may be underestimated because we do not 
take into account the compression of the diffuse gas in the 
spiral structure. Comparing with molecular clouds, we can 
estimate that, having atomic gas concentrated in the spi- 
ral arms would result in a brighter FIR emission, with 
a slighty warmer color temperature since the dust grains 
would be closer (in average) to the heating sources. A pre- 
cise estimate of the magnitude of the effect is beyond the 
scope of this paper. 

The global infrared excess for the model galaxy, IRE, 
is defined as the luminosity ratio IRE= Li2~iooMni/^iyj 
with Llj, = N^yJiVLy and hv^y— 13.6 eV. At the disk 
scale, the IRE takes the value 5.9, in agreement with ob- 



servations of Galactic HII regions (Caux et al. 1985, Myers 
et al. 



1986) 



The diffuse and dense gas (atomic and molecular) have 
the following contributions to the total luminosity of the 
C"*" 158 fim line: 77% from the dense phase and 23% from 
the diffuse phase. Less than 10^ Lq comes from HII re- 
gions. The total emission of the galaxy is 2.5 x 10^ Lq, 
a factor 2.5 lower than the measured value, 6.3 x 10^ 
19931 ). Compared to the 60-100 nm 
the C+ line represents 0.21% of the 
FIR (60-100 fj.m) emission. This figure is comparable to 
the observed ratio for other spiral galaxies with 0.1 - 1 % 
(Lord et al. |1996| ). Nevertheless, the va lue fo r NGC 6946 

and in the 



Lq (Madden et al 
far infrared emission 



199^ , 



was found to be 0.6 % (Madden et al. 
Galaxy, Shibai et al. ( |199lD and Wright et al. ( |199lD have 
measured Lp+ = 0.7 % Lfir with the same definition of 
LpiR as above. 



4-3. Radial profiles 

The 60 /zm radial profile is shown on Fig. 8a. There is a 
large decrease from the inner to the outer parts of the disk, 
about two orders of magnitude. In NGC 694 6, th e same 
behaviour has been observed by Tuffs et al. (199£) using 
ISO. Averaged over the model, the Seo/Sioo infrared color 
appears to be slightly different in the two phases: 0.32 
for the diffuse phase and 0.38 for the dense phase. This 
FIR color decreases with increasing radius from 0.40 in 



the center to 0.23 at R ~5 kpc (Fig 8b), in agreement 
with the maps by Engargiola (1991). The decrease is seen 
in both phases, with Sgo/Sioo ranging from 0.35 to 0.23 
for the diffuse phase, and from 0.42 to 0.28 for the dense 
phase. 

The radial profile of the intensity of the C'^ ^^3/2 - 
158/^TO line (Fig. 8c) shows a much flatter gradient 
than the FIR emission. This is due to the logarithmic de- 
pendence of the line intensity on the incident radiation 
field in PDRs. As shown in Fig. 8c and 8d, the diffuse 
atomic gas is the main source of C"^ emission at large dis- 
tance from the nucleus, for radii larger than 4 kpc. It is 
thus possible to determine the intrinsic L(^+ /L pm lumi- 
nosity ratio from the two gas phases, using the data at 
R ~ 2kpc for the molecular gas and data at i? > 5kpc for 
the atomic gas. We find that Lc+fLpjn is equal to 0.10 
% in the diffuse phase and to 0.25% in the dense phase. 



4.4. Maps 

We show on Fig. 9 face-on maps of 100 /im, C+ and 
UV(912-2000 A) emissions. Edge-on maps at the same 
wavelengths are shown in Fig. 10 for comparison. Com- 
pared to the C+ ob servat ions of the edge-on galaxy NGC 
891 (Madden et al. 1994), there is an overall agreement. 
In particular, the scale height in 0+ is predicted to be 
larger than the scale height of the CO emission, due to 
the contribution of the diffuse neutral and ionized media 
which have a larger scale height (Fig. 11). 

In the face-on C'^ map, there is a large hole in the inter- 
arm regions in the NW, at a similar position to the hole de- 



tected by Madden et al. (1993) with the KAO. This hole is 



due to the lower density of molecular gas and of OB associ- 
ations in the interarm regions. Therefore few UV photons 
illuminate this region and the radiation field is very low. 
The map shows many details and a large contrast between 
arm and interarm regions. We have smoothed the image 
from the model to the resolution of the KAO observations 
(50" beam = 1.2 kpc at the distance of NGC 6946). The 
contrast between the brightest regions and the disk drops 
by a large factor. This resolution effect may explain the 
low dynamical range found in the observed data. If PDRs 
are the main source of C"*" 158/im radiation in galaxies, 
we predict that the emission should have more contrast 
at higher spatial resolution. This could be tested by maps 
of external galaxies made with the future Stratospheric 
Observatory For Infrared Astronomy (SOFIA). 

The edge-on maps at 100 /im and in C^ are fairly 
symmetrical with respect to the center. The edge-on C"*" 
map shows however a hole in the central region ( r < 
500pc) which does not appear on the 100 /xm map. This 
hole is largely due to the large opacity for these lines of 
sight (rc+ = 0.4). 
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4-5. Sensitivity of the model to input parameters 

The model results are of course sensitive to the input pa- 
rameters, therefore we have run different models deviating 
from the standard model by one parameter. 

Because of the poor knowledge of the albedo in UV, we 
have run a model with a lower albedo of dust grains, u}= 
0.4. Wo find that the opacity increases to 1.0 at 1000 A & 
0.90 at 2000 A. The 60 l^m emission from the dense phase 
increases by 5%, while the 100 and 200 /xm emissions both 
decrease by 10%. This difference in far infrared emission 
is due to the moderate increase of the opacity which leads 
to a warmer dust temperature. The effect on the emission 
from the diffuse phase is negligible. 

A more extreme case is for a null albedo, suppressing 
any scattering effect. In that case, wc maximize the UV 
opacity and the FIR emission. The opacity increases to 
1.25 at 1000 A and 1.01 at 2000 A respectively. As a rcsuh 
of this larger absorption, the 60-200/xm emission increases 
by 47 % as compared to the standard model. 

In another run, wc have kept the total mass of molec- 
ular gas constant, but used a lower mean density, 20 H2 
cm~^ instead of 50 H2 cm~^, to increase the clouds sizes. 
The volume filling factor is then 3.8 %. These larger clouds 
block more light, and 30 % only of the molecular cells are 
heated, instead of 40% in the standard model. As a result, 
the 60-200 fj,m luminosity decreases by 10%, to 16.7 x 10^ 

We have also investigated the effect of the number of 
OB associations: we have kept the same star formation 
rate but have gathered adjacent associations to form more 
powerful sources. As a consequence, hqb decreases from 
12000 to 3000. Then a smaller fraction of the cloud popu- 
lation is heated, 15%, as compared to 40% in the standard 
model. But because these cells are heated by more power- 
ful OB associations, the far-infrared emission is larger and 
reaches 20.4 x 10^ Lq. Thus the FIR emission depends 
slightly on the number of associations. The C+ emission 
decreases to 1.8 x 10^ Lq, because of the smaller number 
of illuminated clouds. 

If we now increase the SF efficiency, from 5 to 10%, 
so as to double the UV luminosity, the production rate 
of Lyman continuum photons increases by 80 %. In that 
case, the mean UV opacity is 0.78. The FIR 60-200 fim 
luminosity increases by 55 % to 28.7 x 10^ Lq. This shows 
that the FIR emission is not a linear function of the UV 
luminosity in our model. This non-linear behaviour arises 
because the opacity is largely controlled by geometrical 
effects. With a larger star formation activity, HII regions 
are very large and can destroy molecular clouds efficiently. 
Thus the mass of molecular gas decreases in the model 
with a higher SFR. This is the main reason for the non- 
linear behaviour. This result has been established with the 
same number of OB associations, while an increased SFR 
will probably lead to more associations in the disk. How- 



ever we have previously shown that the FIR emission does 
not depend strongly on the number of OB associations. 

We have investigated the effect of the atomic density 
on the size of HII regions, because we probably overes- 
timate the diameter of HII regions, using a mean atomic 
density and neglecting the dust absorption. If the local gas 
density is multiplied by two, the volume of the Stromgren 
sphere is 4 times smaller than in the standard model. The 
60-200 fim luminosity of dense molecular gas increases by 
10% to 1.1 X 10^° Lq. This is explained by the reduced 
destructive effect of HII regions on molecular clouds, and 
then the larger chance for photons to be absorbed by 
molecular gas. The respective contributions from the dif- 
fuse and dense gas to the FIR(60-200 /um) are now 34% 
and 66%. 

This last test shows that the distance between clouds 

and OB associations has a strong influence on the UV re- 
processing by dense gas. For the standard model, we have 
calculated the mean distance between an OB association 
and the nearest cloud edge, dps /cloud, slucI have found a 
value of 35 pc, the mean distance between clouds centers 
is 37 pc. To have a larger separation between clouds and 
OB associations, wc have increased ^escape to 30 kms"-"^. 
We obtain doB/cioud= 39 pc. The FIR(60-200 jim) emis- 
sion from the dense phase decreases by 15% because of 
the smaller solid angles of the clouds viewed from the as- 
sociations. As for the FIR (60-200 jim) from the diffuse 
phase, it slightly increases by 6%. 

We have shown that part of the UV opacity is due to 
geometrical effects. Indeed the UV opacity is lower when 
the molecular clouds are distributed uniformly in the disk. 
We have used an earlier epoch of the simulation, when the 
distribution of gas clouds is axisymmetric. We have kept 
the same value for the other parameters (number of OB 
associations, star formation rate, etc.). In that case, the 
clouds occupy a larger fraction volume of the galactic disk, 
and the mean distance between clouds increases. Because 
of this larger mean distance between clouds, the opacity 
at 1000 A decreases to 0.47. 

5. Discussion 8z conclusions 

We have shown that with simple assumptions about the 
birth of massive stars and their relationships with the 
ISM, we can reproduce qualitatively and quantitatively 
the characteristics of the UV, Ha and FIR emissions of 
a particular object, the Sc galaxy NGC 6946. For such a 
galaxy with a prominent spiral structure, having a large 
mass of neutral gas, and forming stars actively, the ob- 
served far infrared emission is produced both in molecular 
gas and in the diffuse atomic gas. More precisely, 54 % of 
the FIR (60-200 /xm) emission comes from dust grains in 
giant molecular cloTids. Dust in the diffuse ncTitral atomic 
gas contributes to about 46 % of the total FIR luminosity. 

We have evaluated the respective contributions of the 
UV radiation from massive stars and of the radiation field 
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from the old stellar population. Wc find that 72 % of 
the FIR luminosity can be attributed to UV heated dust 
grains, which reside mostly in molecular clouds envelopes. 
The remaining 28 % is due to dust heated by the radiation 
field from the old stellar population at locations far away 
from OB associations. 

We have calculated the emission of the model galaxy 
in the ^^3/2 — ^ P1/2 fine structure line of 0+ at 158 (xm. 
In the spiral arms, photon dissociation regions at the sur- 
faces of molecular clouds are the main source of the emis- 
sion. We have found a large arm-interarm contrast in this 
line. This effect could be tested by high angular resolution 
maps of galaxies. It results naturally from the combina- 
tion of a lower gas density and lower radiation field in the 
interarm regions, because of the short mean free path for 
UV photons, '-^ 440 pc. As a whole, PDRs represent 76% 
of the emission. The contribution from the diffuse phase 
is found to be ~ 24 %. Our model is able to account for 
about 40 % of the observed 0+ emission of NGC 6946. 
The emission from PDRs should be viewed as a lower limit 



since we use the model by Le Bourlot et al. ( 1993 ) with 
low abundances of carbon and other elements in the gas 
phase: [C]/[H]= 3 x 10^*^^ The average val ue is 1.3 x lO"'' 
for Galactic diffuse clouds (Snow & Witt 1996| ), a factor 
of 4 larger than the value used in the model. Since the C+ 
158 /im emission scales roughly with the column density, 
hence the carbon abundance, the total C"'' luminosity from 
PDRs could be larger by at least a factor three than our 
current model prediction. This would increase the contri- 
bution from PDRs to the total C+ emission of the model 
galaxy: with this scaling factor, the predicted C+ lumi- 
nosity of PDRs would reach 6 x lO^L©. Moreover, the C"*" 
emission from the diffuse gas is overestimated, because the 
model of Wolfire et al. ( |1995D assumes [C]/[H] - 3 x 10"^ 
in the gas. Thus the diffuse emission could be 2-3 times 
smaller than in our standard model. Adopting [C]/[H] = 1.3 
X 10^^ in both phases would thus enhance the differences 
of Lqjt / Lpjft between dense and diffuse gas. 

The knowledge of the cloudy nature of the ISM, and 
of the global structure of the galaxy, is important to de- 
termine how far UV photons can travel away from OB 
associations. The filling factor and the mass/radius scal- 
ing law appear to be major parameters for the transfer 
of stellar radiation in the galaxy disk, because they de- 
termine at the same time the obscuration and the size of 
the emitting regions. Other important parameters are the 
number of OB associations and the sizes of HII regions, 
because with a large number of OB associations or with 
small HII regions, molecular clouds are on average closer 
to massive stars, and are thus more efficiently heated. 

In all the models we ran, we have found that the aver- 
age internal UV opacity is of the order 0.8. The discs are 
therefore mo derat ely opaque in the UV, as measured by 
Buat & Xu (|1996| ). This moderate opacity holds for face- 
on discs. Edge-on discs are quite opaque, with a small 
fraction of the Imninosity escaping, less than 1.0% of the 



face-on luminosity. This fraction corresponds to an equiv- 
alent extinction of 5 magnitudes in the UV. 

These results have been obtained using a crude de- 
scription of the interstellar medium. The adopted spatial 
resolution results from a compromise between astrophysi- 
cal requirements and computational needs, but is certainly 
very poor compared to the complexity of the interstellar 
medium. The good agreement of the observed and pre- 
dicted large scale properties shows nevertheless that the 
transfer of UV radiation, and the role of the radiation for 
the gas and dust heating, are correctly described at the 12 
pc scale. This is in agreement with previous works estimat- 
ing that dust heating by UV radiation occurs principally 



at large distances fr om m assive stars (Murthy et al. 1992 
Leisawitz & Hauser 1988). 
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Figure captions 

Fig. 1. C0(J=2-1) map at 13" resolution obtained with 
the IRAM 30ni radiotelescope. 

Fig. 2. V-I image, obtained with the 1.20m telescope at 
the Observatoire de Haute Provenee (P. Boisse, private 
communication). The field of view is 10' by 10', the pixel 
size is 2.3" . The gray scale runs from white for blue colors 
to black is for red colors. 

Fig. 3a. R band image of NGC 6946 used as input for the 

calculation of the gravitational potential. Field stars have 
been removed, the image has then been rotated and de- 
projected to get a face-on view of NGC 6946. 
Fig. 3b. A face-on view of the model galaxy, with molec- 
ular clouds drawn as circles and OB associations drawn 
as stars. Only 20 % of the OB associations and molecular 
clouds are drawn. 

Fig. 3c. A close-up view of Fig. 3b. The clouds are drawn 
at their exact size in the model. 

Fig. 4. Adopted rotation curve (dot-dashed line), angular 
frequency fl{r) (full line), and ± f curves. The corota- 
tion for the adopted pattern speed is located at r = 3.5 
kpc. 

Fig. 5. Distribution of the production rate of Lyman con- 
tinuum photons NLyc, for the modelled population of OB 

associations (full line) in NGC 6946, and for the Galaxy 
(dot-dashed curve) (McKee & Wilhams (1997)). The plot 
can be read for instance as 500 OB associations out of 
12000 have N^yc > 10'™ s-^ 

Fig. 6. Distribution of the distance from their parent OB 
association, travelled by UV photons before absorption. 
The median distance is 120 pc, and the mean distance is 
440 pc. 

Fig. 7. Distribution of the UV intensity measured rela- 
tive to the ISRF in the modelled galaxy, for cells filled 
with diffuse gas (dot-dashed line) and for cells filled with 
molecular gas (full line). 

Fig. 8a. Radial profiles of the 60 /im surface brightness. 
The thin solid line shows the combination of the two gas 
phases, while the triangles show the contribution of the 
diffuse atomic gas, and the squares the contribution from 
the molecular gas. The thick solid line shows the mean 
radial profile deduced by Engargiola (1991) from an IRAS 
image. 

Fig. 8b. Radial profiles for the ratio of 60 and 100 /xm 
fluxes, global and for the two phases. 
Fig. 8c. Radial profiles of the C+ emission: global and for 
the two phases. 

Fig. 8d. Radial profiles of L^/-'^(60-i00(um)- global and for 

the two phases. 

Fig. 9. Face-on views of: 

a) 100 fim emission, at 48 pc resolution. The gray scale 
ranges from 1 to 10^ Lq pc~^. We have overlaid contours 
of the same map convolved with a 750 pc beam: the levels 
are at 10, 30, 60, ICQ, 200 Lq pc'^ 

b) C+ line, at 48 pc resolution. The gray scale ranges from 



10~^ to 1 Lq pc~^. Overlaid contour levels from 0.1 to 0.6 
by 0.1 L0 pc~^ for the same image convolved with a 750 
pc beam. 

c) emergent UV surface brightness (resolution: 48 pc). The 
gray scale ranges from 10~^ to 10"'' Lq pc~^. 
Fig. 10. Edge-on views of NGC 6946, at 48 pc resolution. 
The linear scale is not identical for both axes. 

a) 100 /iTO IRAS band, with contour levels at 10^, 3 x 10^, 
6 X 10^ and 10'' Lq pc'^. 

b) C+ , at the same resolution, accounting for the line opac- 
ity. The gray scale ranges from to 50 Lq pc~^. Contour 
levels from 10 to 50 by 10 Lq pc^^. 

Fig. 1 1 . Average vertical profiles through the disk for the 
C+ (dot-dashed line) and C0(l-O)(full line) emissions. 
The CO(l-O) profile presents a smaller scaleheight (40 pc) 
than the C+ profile (100 pc). We have assumed that the 
CO(l-O) emission is proportional to the molecular gas col- 
umn density. 
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Table 1. Observed parameters for NGC 6946. 



Observed data 


Value 


Distance 




5 Mpc 


Inclination 




34° 


Mh2 (R < 6') 


2.0 


10" Mq ^ 


Mhi (R < 6') 


2.6 


10" Mq 2 


Ll950-2050A 


5.4 


10» Lq 3 




1.0 


10» Lq * 


B? 


8 


49 mag 


IRAS 12^m 


2.3 lO'* 


Lq ± 20 % ^ 


IRAS 25^m 


1.4 lO'* 


Lq ± 20 % 5 


IRAS GOfim 


4.5 lO'' 


Lq ± 20 % 5 


IRAS lOO^im 


6.6 lO'' 


Lq ± 20 % 5 


IRAS 200^im 


3.4 10" 


Lq ± 20 % 5 


Lc+ 


6.3 


10^ Lq e 


{^thermal 


83 mJy ± 25 ^ 



Reference 



De Vaucouleurs (1979) 
Considere & Athanassoula (1988) 
Young & ScoviUe (1982) 
Boulanger & Viall efond (1992) 
Buat et al. ( |l989|) 
DeGioa et al. ( |1984| )) 
RC3 

Engargiola (1991) 
Engargiola (1991) 
Engargiola (1991) 
Engargiola (1991) 
Engargiola (1 991) 
Madden et al. (|1993| ) 
Klein et al. ([I982|) 



We use the value 3.8 x 10' 



'^^ erg s ^ for the solar luminosity at any wavelength. 



We use N//2/I 



CO{1-0) 



2.3 10^° mol cm /(K kms ^) and a total intensity I 



CO{1-0) = 



569 Kkms in a 45" beam, we do 



not account for projection effects. 

^ We use Mhi = 1.9 x 10^° Mq at a distance of 10 Mpc, corresponding to 4.7 x lO" Mq for a distance of 5 Mpc. By performing 
an integration over the HI map, we estimate that the part with R < 6' contributes to 55 % of the global emission. 
^ The 2000 A flux, 6.92 x 10"^^ erg cm ^s ^, is corrected for the Galactic extinction and integrated over a 100 A band. 



Note that Donas & Deharveng (1987) have reported a flux at 2000 A twenty times lower (1.66 x 10 ^'^ erg cm ^s 



Measurements at other UV wavelengths are available in Rifatto et al. (1995). 
* We use th e uncorrected flux iHa~ 1.31 x 10~^° erg cm~^s~^ and the relationship Lna 



3.13 X 10 

(1996). We have subtracted the contribution from the nucleus, estimated to be 20 % of the total luminosity. Note that 



Dmpc fea from Young 



et al. 

Young et al. ( |199^ ) give a lower total flux, 3.38 x 10 erg cm s ^ , while Kennicutt (|1989|) reports a mean surface brightness 



of 4.2 X lO'^'^ erg s^^ pc"'^ corresponding to a total Ha luminosity of 3 x 10^ Lq at the adopted distance. 
^ We estim ate th e luminosity in the band centered on frequency u, by uFv We use the IRAS and KAO fluxes reported by 
Engargiola (L991) for a radius R < 5.6' and correct from the contribution of the nucleus (central 45") to obtain 12 Jy at 12^m, 
15 Jy at 25/^m, 114 Jy at 60/Ltm, 283 Jy at lOO^m and 288 Jy at 200/^m. The global color for the disk is Seo/Sioo 
far infrared luminosity is Lpjuc 60-ioonm)= Leo^ + Liooa. = 1-1 x 10^° Lq, Lpjjn 60-200^.™) = 1-45 x 10^° Lp 



0.40. The 



We have excluded the contribution of the nucleus and rescaled the dat a by M adden et al. (1993) for the adopted distance. 
^ We use an estimation of 625 mJy for the total flux at 6cm (Klein et al. 1982), subtract the nuclear contribution estimated to 



be 33% of the total flux, and we keep 20% of the resultant disk flux as the thermal component. 



Table 2. Input parameters for the standard model. 



Parameter 


Adopted value 


cell resolution 


12.2 pc 


HI half thickness 


180 pc 


H2 scaleheight 


65 pc 


mean Vos 


10 kms'^ 


n-y 


100 


time step 


lO'' years 


number of molecular clouds 


18000 


number of OB associations 


12000 


star formation rate (2 - 60 Mq) 


4.0 Mq yr"^ 


volume filling factor of H2 


2.5 % 
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Table 3. Output values from the simulation for the standard model. 



Parameter 


Total 


Dense 


Diffuse 


uv 


ISRF 


L912-20OOA (Lq) 


1.3 10"' 


— 


— 


— 


— 


Ll950-2050A (Lq) 


5.2 10** 


— 


— 


— 


— 


L/fa (L0) 


1.2 10** 


— 


— 


— 


— 


B? (mag) 


8.96 


— 


— 


— 


— 


L12 (L0) 


2.1 lO'' 


1.0 10® 


1.0 10® 


1.7 10® 


0.3 10® 


L25 (Lq) 


2.0 10^ 


1.0 10® 


1.0 10® 


1.7 10® 


0.3 10® 


Leo (Lq) 


5.1 10® 


2.9 10® 


2.2 10® 


4.2 10® 


0.9 10® 


Lioo (L©) 


8.5 10® 


4.5 10® 


4.0 10® 


6.1 10® 


2.4 10® 


L200 (Lq) 


4.9 10® 


2.5 10® 


2.5 10® 


3.0 10® 


1.9 10® 


L60-200 (L©) 


18.6 10® 


9.9 10® 


8.6 10® 


13.3 10® 


5.2 10® 


Lc+ (Lq) 


2.6 lO'' 


2.0 lO'' 


0.6 10'' 






Seem (mjy) 


124 










Lc+/Lf/H (60-100) 


0.18 % 


0.25% 


0.10 % 







All luminosities are emergent luminosities. The total UV luminosity generated by 
OB associations in the disk between 912 and 2000 A is 3.10 x 10^° L©. 
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Abstract. We present numerical simulations of radiative 
transfer in the spiral galaxy NGC 6946. The interstel- 
lar medium is represented as a two phases medium, with 
molecular clouds and a smooth diffuse phase. The molec- 
ular gas distribution is calculated in a self-consistent way 
from the distribution of an ensemble of molecular clouds 
evolving in the gravitational potential of NGC 6946. We 
simulate star formation by creating OB associations in 
molecular clouds. The transfer of UV radiation is calcu- 
lated in the clumpy interstellar medium, to determine the 
local UV illumination of molecular clouds. We compute 
the emergent intensity in the UV continuum (912 - 2000 
A), in the Ha and C+ ^P'iji—^ Pi/ 2 lines as well as in the 
continuum at far infrared wavelengths, 60, ICQ & 200 /im. 

It is possible to obtain a consistent picture of this 
galaxy with a global star formation rate of 4 M0yr^^ (for 
stars with masses in the range 2-60 M©) occuring mostly 
in the spiral arms. The close spatial association of massive 
stars and molecular clouds has a profound impact on the 
transfer of UV radiation in the galactic disk and on the 
dust emission. The median distance travelled by UV pho- 
tons is about 120 pc. However, when they have escaped 
from the vicinity of their parent OB associations, UV pho- 
tons may travel quite far in the disk, up to 1 kpc. The UV 
opacity of the model spiral galaxy disk, observed face-on, 
is 0.8 at 1000 A and 0.7 at 2000 A. 

For radii less than 4 kpc, the C"*" 158 /xm line is mostly 
produced in photodissociation regions at the surfaces of 
molecular clouds. The C+ emission from diffuse atomic 
gas accounts for about 20% of the total. It becomes sig- 
nificant at large distance from the nucleus (r > 4 kpc). 
Molecular clouds and diffuse atomic gas have almost equal 
contributions to the total far infrared emission from 60 to 
200 iim. As a whole, 72% of the 60 - 200 //to FIR emis- 
sion can be attributed to dust grains heated by the UV 
radiation of massive stars and 28% by the radiation field 
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of the old stellar population. 



Key words: Galaxies : ISM, Galaxies : individual: NGC 
6946, - ISM : molecules, ISM : dust 



1. Introduction 

Because the peak of the thermal dust emission occurs 
in the far-infrared (60-200 /im) for Galactic interstellar 
clouds, this wavelength domain is recognized as a very im- 
portant 'window' to study the interstellar medium, both 
in the Milky Way and in external galaxies. However, due 
to the lack of spatial resolution provided by telescopes 
operating at these wavelengths, the interpretation of the 
signals measured in external galaxies is still a subject of 
controversy. The debate is focussed on the origins of the 
heating of dust grains and of the "^Pz/-! P1/2 C"*" line at 
158 /tm. This fine structure line of ionized carbon is one of 
the main cooling lines of neutral atomic gas (Wolfire et al. 
1995), but it is also seen in HII regions, with different exci- 
tation mechanisms in these two cases. Other sources of C+ 
emission are the dense photodissociation regions (PDRs) 
where intense UV radiation from massive stars impinges 
on molecular clouds. In these regions, the photodissocia- 
tion of carbon monoxide creates a layer of ionized carbon 
at the edge of the molecular cloud, which is a source of 
intense C"*" emission. 

The far infrared emission is produced by dust grains 
heated by star fight but the relative contribution of young 
massive stars on one hand, and of the bulk of the stel- 
lar population on the other hand, are still debated (see 
e.g. Walterbos & Greenawalt 1996, Devereux & Young 
1990, Cox & Mezger 1989, Thronson et al. 1990, Rice et 
al. 1990). Using low resolution data from IRAS, IR col- 
ors have been used to assess the respective roles of star 
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forming regions and cirrus in external galaxies at a global 
scale (e.g. Calzetti et al. 1995). It is difficult to find spi- 
ral galaxies dominated by a single population of heating 
sources, since the relative contributions of massive stars 
and of the disk population vary from galaxy to galaxy, 
and probably from place to place in a given galaxy. A 
main problem of these studies is the lack of spatial resolu- 
tion at far infrared wavelengths. Even when it is possible 
to use high spatial resolution data at millimeter and sub- 
miUimeter wavelengths, the gap in spatial resolution from 
the local interstellar clouds to external galaxies is huge: 
for nearby galaxies located at distances £> of a few Mpc, 
a one arcmin beam encompasses 0.3(r'/lMpc)kpc, much 
larger than the size of a single molecular complex in the 
Galaxy, whereas structures down to 0.01 pc are commonly 
observed in local clouds (Falgarone et al. 1992) and in the 
Magellanic Clouds (Rubio et al. 1993). 

Numerical simulations provide a way out of this prob- 
lem: it is possible to reproduce the observed molecular 
gas distribution of nearby spiral galaxies with numerical 
simulations using observed data as input parameters. For 
example Garcia-Burillo et al. (1993) were able to fit the 
spatial distribution and kinematics of the molecular gas in 
M 51 using the cloud colhsion code developped by Combes 
& Gerin (1985). Though this code has a spatial resolution 
of a few hundred parsecs (the cell size for the large scale 
dynamics), it is possible to include "micro-physics" at the 
parsec scale inside each cell. We have taken this approach 
and implemented star formation in this code to study the 
far infrared and C"*" emissions in the spiral galaxy NGC 
6946. The next section summarizes the current data on 
NGC 6946. Wc describe the model in Sect. 3 and present 
the results for NGC 6946 in Sect. 4. The imphcations of 
this work for the interpretation of the C+ and FIR emis- 
sion of spiral galaxies are discussed in Sect. 5. 

2. NGC 6946 

NGC 6946 is a nearby Scd galaxy with a low inchna- 
tion, i = 34° (Considere & Athanassoula 1988). At the 
adopted distance of 5 Mpc, 1 arcmin on the sky corre- 
sponds to a linear size of 1.5 kpc. The main characteris- 
tics of NGC 6946 are summarized in Table 1. It shows a 
weak bar (Martin 1995) and a prominent spiral structure 
with both an m = 2 and an m = 4 pattern (Considere & 
Athanassoula 1988). With this open spiral structure, the 
arm and interarm regions are resolved at moderate spatial 
resolution (30"). The gas distribution is well known since 
NGC 6946 has been mapped at high spatial resolution in 
HI (Boulanger & Viallefond 1992), and CO (Casoh et al. 
1990, Clausset et al. 1991). Figure 1 presents a C0(2-l) 
map obtained at the IRAM 30m telescope (Clausset et al. 
1991), and Fig. 2 a V-I color map obtained at the Ob- 
servatoire de Haute-Provence (OHP) (P. Boisse, private 
communication). NGC 6946 is marginally resolved in the 
far infrared maps obtained by IRAS and the KAO (Engar- 



giola 1991). It also belongs to the normal galaxy sample 
studied hy ISO and as such has been extensively mapped 
in the mid and far infrared (Malhotra et al. 1996, Helou 
et al. 1996, Lu et al. 1996, Tuffs et al. 1996). NGC 6946 is 
forming stars actively, as revealed by its bright Ha emis- 
sion (Bonnarel et al. 1988a, 1988b, Kennicutt 1989) and 
the numerous supernovae (Li & Li 1995). Madden et al. 
(1993) have mapped the C"*" 158 /xm emission at 45" reso- 
lution with the KAO. Bright C"*" emission is found in the 
nucleus and in the disk. At this scale, the C+ emission is 
correlated with the CO emission in the spiral arms and 
with the HI emission in the outer disk. 

3. A model of NGC 6946 

3.1. Representation of the interstellar medium. 

We have chosen to include two phases for the interstellar 
medium : 

i) a dense molecular phase composed of spherical molec- 
ular clouds with masses ranging from 10'^ to 10® M0 and 
mean density 50 H2 cm~^. This mean density is used to 
determine cloud sizes from their masses. The half power 
width of the molecular layer, defined as < z'^ >^/^, is equal 
to 65 pc. The global volume filling factor of this molec- 
ular phase in the disk, (H2), amounts to 2.5% . The 
clumpy structure of molecular clouds is taken into account 
by using a larger density, 5 x 10^ H2 cm~^, to compute the 
C+ emission. The C"*" emission is not sensitive to the gas 
density when it is larger than ~ 10^ cm^^ (see Tielens 
& Hollenbach 1985 and Sect. 3.6). Indeed, 0+ observa- 
tions of nearby molecular clouds show that C+ emission 
is detected over the whole projected surface of molecular 
clouds (Jaffe et al. 1994). The total mass of dense molecu- 
lar gas in the model is set to match the value deduced from 
CO(l-O) observations with a conversion factor of 2.3 x 10^° 
H2cm-2 (K kms-i)-i, 2.0 x lO^^ Mg (Table 1). 

ii) a neutral diffuse phase, with a nearly constant HI den- 
sity. We have adjusted the value of the mean HI density us- 
ing the data from Boulanger & Viallefond (1992): it varies 
slowly from 0.8 H cm^^ at the center to 2 H cm^^ in the 
outer disk (r ~ 6 kpc) . We have assumed that the HI disk 
has a vertical thickness of 360 pc = 2iJ^, and that the 
density is uniform from — i?^ to The velocity dis- 
persion of the diffuse gas is set to 10 kms~^. The total 
mass of diffuse neutral gas inside a radius of 6 arcmin is 
taken from the HI observations to be 2.6 x 10^ M©. 

In addition to these two phases, we take into account 
the formation of HII regions around OB associations (see 
below). For the numerical calculations, we represent the 
galaxy disk as a large 3D grid. The total grid length L 
and the cell size rceii are chosen to match the following 
constraints: 

- the grid size has to be at least as large as the optical 
disk, which has a diameter of 12.5 kpc, in order to be able 
to follow the path of photons in the external regions of the 



S. Sauty, M. Gerin , F. Casoli: FIR and C*" emissions of spiral galaxies disks. 



3 



galaxy, and the heating of the neutral gas at large distance 
from the nucleus. 

- the spatial resolution, rceii, should be good enough to 
describe the clumpy structure of the medium: in partic- 
ular the smallest clouds in the ensemble should occupy 
at least one cell. With the adopted parameters for the 
clouds, the radius of the smallest clouds with a mass of 
10^ M0 is 4.5 pc. This number is only an estimate of the 
actual size of a given cloud. It is possible to use an alter- 
native method based on the scaling relations for molecular 
clouds, the so-called Larson's laws. Using the mass-radius 
relation M(Mq) = 100R'^{pc^) (Falgaroneet al. 1992), the 
radius of a 10^ Mq cloud is then 3 pc. Another constraint 
on the spatial resolution is provided by the size of OB as- 
sociations, which is generally larger than 30 pc (Garmany 
1994) and their distance from molecular clouds. Leisawitz 
(1991) has estimated the mean distance between OB as- 
sociations and their parent cloud to be about 50 parsecs. 

We have chosen as the best compromise to represent 
the galaxy with a 1024 x 1024 x 24 cell structure. The 
galaxy size is then L = 2 x Rmax = 12.5 kpc and the 
resolution rceii = 12.2 pc. The cells are filled with neu- 
tral atomic gas, or with molecular gas at the molecular 
cloud positions, an 10^ M0 molecular cloud then occu- 
pies 3x3x3 cells. Finally, we include the ionized gas 
in the Stromgren spheres centered on each OB associa- 
tion, and replace the molecular and atomic gas with ion- 
ized gas whenever required. The radii of the Stromgren 
spheres are calculated assuming classical ionization in HI- 
bounded HII regions (Miller & Cox 1993). Because of the 
coarse spatial resolution of the dynamical code (200 pc), 
we do not attempt to reproduce the interstellar medium 
in the central region of NGC 6946 (r < 500 pc). 

3.2. Spatial distribution of molecular clouds 

This section of the code uses the cloud-cloud collision code 
described by Combes & Gerin (1985), and used by Garcia- 
Burillo et al. (1993) and Gerin et al. (1991) to model the 
gas dynamics in nearby galaxies. The molecular clouds 
move in the gravitational potential of the galaxy, they 
grow through cloud-cloud collisions, which are sticky pro- 
cesses, and are disrupted by simulated SN events. The 
molecular gas is immediately recycled into small molec- 
ular clouds. The gravitational potential is deduced from 
an R band image of NGC 6946 (Viallefond & Bonnarel, 
private communication, Bonnarel et al. 1988), with fore- 
ground stars removed. The image has been rotated to put 
the major axis vertical, and deprojected to face-on, using 
as projection parameters PA = 69° and i = 34 °. The 
image has then been binned to 256 x 256 pixels. This R 
band image covers 8.5' x 8.5' on the sky, corresponding 
to a radius of 6.24 kpc at the assumed distance of NGC 
6946, 5 Mpc. The actual spatial resolution amounts to a 
few times the cell size, about 200 pc. The following step 
is to build a good gravitational potential from this image. 



As in Garcia-Burillo et al. (1993), the calculation is done 
in two steps. The axisymmetric part of the potential is 
obtained by assuming a constant mass to light ratio, and 
adjusting this constant to reproduce the observed CO and 
HI rotation curves. The perturbations due to the bar and 
spiral arms are included in the non axisymmetric part of 
the potential. The last parameter to adjust is fip, the pat- 
tern speed of the density wave (bar + spiral arms) . The gas 
distribution and velocity field are very sensitive to ftp. We 
find that ilp= 42 kms~^ kpc~^ gives the best results. The 
corotation is found at a radius of 3.5 kpc, close to the end 
of the bar, and the OLR hes outside of the disk. Figure 3a 
presents the deprojected R-band image. Figures 3b and 3c 
present an example of the obtained molecular cloud distri- 
bution. The gas clouds are concentrated in the spiral arms, 
with few molecular clouds at radii larger than 4 kpc. The 
model rotation curve is shown in Fig. 4, together with the 
angular frequency f2(r). We have not attempted to model 
the compression of diffuse gas in the spiral structure, and 
have kept an axisymmetric distribution of diffuse atomic 
gas. 

3.3. Stellar population 

Once the gas is distributed over the galactic disk, we cre- 
ate OB associations, with star masses ranging from 10 to 
60 M0 (stellar types from B2 to 05). We have chosen not 
to include stars less massive than 10 Mq because their life- 
time becomes a significant fraction of the rotation period. 
Shortward of 2000 A, the UV radiation is mostly produced 
by OB stars and the contribution by late B type stars is 
at most 20% (Walterbos & Greenawalt 1996, Mathis et al. 
1983). The choice of the upper mass cutoff at 60 Mq is mo- 
tivated by the work of Heydari-Malayeri & Beuzit (1994) 
who have shown that suspected very massive stars (~ 100 
Mq) are actually clusters of less massive stars. Also, a 
very short time step is required to sample adequately the 
lifetime of these very massive stars. 

We allow only clouds more massive than 4 x 10" M© 
to form massive stars. In the cells satisfying this criterion, 
the OB associations are born at the outer edge of the 
molecular clouds. We have tested different star formation 
laws, so as to match as well as possible the Ha radial pro- 
files. The best fits are obtained with a star formation rate 
(SFR) depending on the angular frequency and the local 
gas surface density as proposed by Wyse & Silk (1989) : 

SFRir) = enir)agas(r) 

with e the star formation efficiency, 0.{r) the angular fre- 
quency at radius r and crgas(r) the HI+H2 surface density 
at position r in the disk. We find that e = 5% provides the 
best match to the Ha radial profile by Kennicutt (1989). 

The stellar mass distribution inside an OB association 
is drawn from the Initial Mass Function (IMF). We have 
chosen an index close to the Salpeter IMF, as suggested 
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by measurements in Galactic and extragalactic OB asso- 
ciations (Massey et al. 1995) : = M-2-3. 

The stellar associations are born with a mean veloc- 
ity relative to their parent cloud chosen from a Gaussian 
of mean value 10 kms^^ and FWHM ~ 10 kms^^. Those 
figures are in good agreement with observed data by Lei- 
sawitz et al. (1989) in their study of the relation of star 
clusters with molecular clouds. With this value, the aver- 
age motion of OB stars relative to their parent cloud is 10 
pc in 10^ years. Finally, stars die after a time equal to their 
Main Sequence lifetimes from Giistcn & Mczgcr (1983). 
Effective temperatures, stellar luminosities and radii, av- 
eraged over the main-sequence lifetime, are taken from 
Cox et al. (1986). Because we are only interested in broad 
band fluxes and colors, we have calculated all stellar fluxes, 
from UV to B bands, in the black-body approximation. 
Lyman continuum radiation production rates, averaged 
on the Main Sequence lifetime, are also taken from Cox 
et al. (1986) and Giisten & Mezger (1983) to have a co- 
herent set of parameters. We use the parametrization as 
a function of the stellar mass, M: 



^offio — f- 



38.3 + 8.16 X logia{M) 



-0.24 X logfoiM) - 0.41 x loglo{M). 

to compute the thermal radio continuum flux. Ha lumi- 
nosity and Stromgrcn sphere diameters. 

We also include the effects of the radiation of massive 
stars on neutral gas, and allow molecular clouds to be par- 
tially eroded and ionized by the radiation of nearby OB 
associations, if they overlap with the Stromgren sphere of 
an OB association. We compute the radius of this sphere, 
Rhii , assuming that all stars are located at the same po- 
sition and have a global production rate of Lyman contin- 
uum photons NLyc, equal to the sum of the contribution 
of the individual stars, and for case B recombination in a 



3_NLyc^y/3 



In 



diffuse medium of density rijj: Rhii = (4 
this formula, as is the hydrogen recombination coefficient 
(Osterbrock 1989). The mean radius for Rhii is 35 pc. 

To constrain the population of massive stars, we calcu- 
late different stellar outputs, namely we perform a detailed 
calculation of the UV radiation field at 912-2000 A, and 
also compute the global U fiux and U radial profile, as 
well as the thermal radio continuum and Ha emission. 

• 6 cm radio-continuum. 
According to Mezger (1972) and Turner & Ho (1994), the 
radio continuum emission of HH regions for case B recom- 
bination, at an electronic temperature Tg of lO'^ K and 
with 45 % of the ionizing photons being converted into 
Ha photons, is directly related to N^^c by: 



Lyc 



D 



_ 1.1 X 10S°(^)(-9^)2. 

This formula does not include any correction for dust ab- 
sorption of the ionizing radiation within the HH region. 



Current estimations are that nearly 50% of the Lyman 
continuum may be absorbed by dust. 
• Ha. 

The relationship between the Ha luminosity and the pro- 
duction rate of Lyman continuum photons, using the same 
assumptions as above, can be deduced from Mezger (1972) 
and Peimbert et al. (1975): 



3.55 X 10 



-46 



The correction for the extinction may amount to about 1 
magnitude at Ha but it is highly uncertain (Mc Kee & 
Williams 1997, van der Hulst et al. 1988). Because of the 
uncertainties involved in this correction, we have preferred 
not to do it. This allows us to check that the energy is 
conserved with a good accuracy in the simulation. 

The formed OB associations have typical production 
rates of Lyman continuum photons, N^yc, in the range 
10^^-10^^ s^^. The cumulative Lyman continuum lumi- 
nosity function of the population is shown in Fig. 5. Its 
shape is similar to the distribution for Galactic HII regions 
(McKee & Williams 1997). There is a small excess around 
Niyc > lO"*® -^ s~^, dominated by the most massive stars 
in the younger associations, and a deficiency of associa- 
tions with Ni^yc > 10^°-^ s^^. As a whole, the match of 
the two distributions is very good. The distribution of the 
intrinsic UV luminosity of OB associations extends over 
four orders of magnitude from 10^ to lO'^ Lq, i.e. from 
small associations with about 20 stars and a total stellar 
mass of 600 M0, up to large associations gathering 200 
OB stars, including a few 50-60 Mq stars, and having a 
total stellar mass of 6000 M©. 



3.4- Dust properties 

We use the average Galactic extinction curve from Fitz- 
patrick & Massa (1988) at UV wavelengths. For visible 
and near infrared wavelengths, we use the work by Seaton 
(1979). This curve is probably valid on a large scale in the 
diffuse medium of NGC 6946 which has a similar metal- 
licity to the Milky Way. We assume a constant metallicity 
and gas to dust ratio in the disk, and use the average value 
for the Milky Way : = 5.8 x lO^i H cm'^ mag-i 

and Av = 3.1 Eb-v (Bohlin et al. 1978). The extinction 
through a cell filled with molecular gas, with the adopted 
spatial resolution, is 2.5 magnitudes. The dust properties 
have been summarized by Bruzual et al. (1988) and Witt 
& Gordon (1996). Apart from the enhanced absorption in 
the 2175 A bump, the dust albedo is fairly constant at UV 
and visible wavelengths at a; ~ 0.55. We include the coher- 
ent scattering of UV light by dust grains. The anisotropy 
is described using the Henyey-Greenstein (1941) function, 
where the anisotropy parameter is defined as g =< cos9 > 
and 9 is the scattering angle. 
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3.5. Radiative transfer 

The observed properties of external galaxies depend on 
tho propagation of tho stollar radiation in tlio interstellar 
medium. To determine the local radiation field, we fol- 
low the propagation of UV photons (912-2000 A) emit- 
ted from the OB stars in the two phases medium. In 
each cell, we compute a local radiation field. We define 
the local radiation field in the UV, Xuv, relative to the 
mean radiation field in the UV at the Solar radius, XO) 
established by Mathis et al. (1983), so that xuv = 
In the Galaxy and at the solar radius, the InterSteflar 
Radiation Field, ISRF in the UV, has been defined as 

Xo = in J = 47r/g^°2T "^^^^ = ^"^^ ^ cm-^g-i 
and through visible/IR bands as Go = 4tt J^'^'p^.J\dX = 

2 X 10~^ erg cm~^s~^, where is the specific intensity 
of the radiation field averaged over Att sr (Mathis et al. 
1983). 

On their path, photons can be scattered and/or ab- 
sorbed by dust in both neutral phases. We have chosen 
not to follow individual photons which would have been 
time consuming, but to use instead pseudo-photons rep- 
resenting a collection of N photons. We are then able to 
probe a larger area with a lower number of photons. We 
sample the 912-2000 A interval with 20 bins of constant 
wavelength width AA, and launch n-y pseudo-photons per 
association per wavelength bin. The longward limit has 
been set to 2000 A, to avoid having to take into acount 
the contribution from star types later than B to the ISRF 
(Walterbos & Grccnawalt, 1996). 

The pseudo-photons emitted by the OB association 
number i carry a fraction dX of the luminosity L\ radi- 
ated by this OB association. They travel from the center 
of a cell (size rceii) to an adjacent one, and the energy 
absorbed by the interstellar medium, when non-zero, is 
left on the common face of these 2 cells. The absorbed en- 
ergy is reprocessed in the far-infrared. The pseudo-photon 
energy along its path from association number i can be 
written as : 



In this formula, a' = a-' (a;, y, z) is the product of 
the probabilities for non absorption in each cell along the 
travel from the OB association number i to the cell po- 
sition {x,y,z), or with lo{X) the dust albedo and t\ the 
total cell opacity including scattering : 

in the diffuse phase, and 

= Lj{X) 
in the H2 phase. 

When leaving an OB association, the direction for each 
pseudo-photon is uniformly chosen on the unit sphere. 
These pseudo-photons travel in the two-phase interstellar 
medium. There are different possibilities when reaching a 
new cell : 



i) The gas in the cell is diffuse and atomic. Then the 
pseudo-photon can either : 

- be absorbed partially in the cell. A fraction of the lu- 
minosity is left, and the pseudo-photon continues in the 
same direction with a lower luminosity. 

- be scattered. 

ii) The cell is filled with molecular gas. 

The pseudo-photon is partly absorbed, and partly back- 
ward scattered, the ratio between the energy left in the cell 
and the total energy of the pseudo-photon depends on the 
albedo as 1 — w(A). We have chosen backward scattering in 
that case to prevent the pseudo-photon from interacting 
with the same molecular cloud several times. 
After a scattering event, a scattering angle 6 is chosen 
according to the anisotropy function, and the azimuthal 
angle (p is uniformly chosen in the interval [— 7r,7r]. Then 
the new direction is easily deduced from the previous one 
(Witt 1977). 

Since we cannot store all the pseudo-photons' incident 
directions, and in order to define an isotropic radiation 
field in each cell, we assume that the equivalent surface 
of the cell is 6r^g„. Then, the UV radiation field in the 
cell {x, y, z) due to the contributions of all pseudo-photons 
travelling through this cell, can be expressed relative to 
the Galactic ISRF xo '■ 

Because of the coarse angular resolution and the low num- 
ber of pseudo-photons leaving each OB association, the 
resultant map of the radiation field presents strong fiuc- 
tuations. We have chosen to smooth the map of the ra- 
diation field by averaging the data in nearby cells, typi- 
cally 3x3 cells. Furthermore, some cells are never visited 
by UV photons, for example in the interarm region or in 
the outer disk (r > 5 kpc). In that case, we use as incident 
radiation field, the Galactic ISRF (longward of 2000 A), 
scaled by the local surface brightness in the R band image 
to take into account the radial variation of the radiation 
field from the old stellar population between the central 
regions (r < 1 kpc) and the outer disk (r > 5 kpc). The 
reference value is given in units of 47r f ^^™o J.dX = Gn - 

•'200071 

Xo = 1-82 X 10~^ erg cm~^s~^. The position for the ref- 
erence value has been chosen at the edge of a spiral arm, 
at a distance r = 4 kpc from the center, where the UV 
radiation field xuv is close to 1. 

3.6. Emergent emission. 

Once the local UV energy density has been calculated, 
models are used to determine the C"*" and FIR emissions 
from each cell, assuming that they arise from the same 
area as the one used for the calculation of the UV radiation 
field. For the dust emission, we use the model by Desert et 
al. (1990), which has 3 different components: PAHs which 
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are fully ionized when xuv > ^, very small grains and big 
grains. We calculate the dust omission in the four IRAS 
bands at 12, 25, 60 and 100 H'ln plus an additional band 
at 200 /xm, as the reprocessing of the combination of the 
UV radiation field (912 - 2000 A) described by xuv and 
of the ISRF for the 2000 A-2 /im part of the spectrum. 
We decrease the PAH abundance in large radiation field 
environments, as suggested by Ryter et al. (1987), to one 
fifth of the standard value when xuv > 100. 

Very large molecular clouds, with masses larger than 
10^ Mq occupy more than 4 cells in the grid. The inner 
cell is not directly exposed to the UV radiation and for 
this cell we assume that the ISRF is attenuated by 2.5 
magnitudes of visual extinction. 

We assume low optical depth in the mid and far 
infrared. From the analysis of the COBE maps of the 
Galaxy, Boulanger et al. (1996) deduce tx/'^h= 1 x 10~^^ 
cm^ (A /250 /im)~^, which combined with the mean 
column density of individual clouds, Njf = 1.5 x 10 
cm^^, gives an opacity of 2 x 10^^ at 200 fim. 

According to the dominant phase in a given cell of the 
model, the emergent infrared emission can be : 

- reprocessing of the whole incident stellar radiation for 
molecular clouds which are totally optically thick in the 
UV. All impinging radiation is completely reprocessed in 
infrared emission from the outer cells of molecular clouds, 

- proportional to the gas column density for the diffuse 
medium, which is optically thin in the UV. We assume no 
UV extinction at the 12 pc scale. 

We use the PhotoDissociation Region (PDR) model by 
Le Bourlot et al. (1993) to estimate the ^Pg/j P1/2 C+ 
emission at the surfaces of molecular clouds. We use a con- 
stant molecular hydrogen density of 5 x 10^ H2 cm^^ and 
the incident UV field xuv ■ The density is not a critical pa- 
rameter as long as it is higher than the critical density for 
coUisional excitation of C"*" (1000 H cm~^) (see Tielens & 
Hollenbach 1985). Furthermore, we assume that the whole 
surface of clouds contributes to the C"*" emission. As for 
the C"*" emission from the diffuse neutral phase, we use 
the model by Wolfire et al. (1995) for a two phase neutral 
atomic medium. Finally, the contribution from the ionized 
gas in the HII regions around the OB associations is also 
included. We assume that the gas has the same density as 
the diffuse medium, rig = uh and an electronic tempera- 
ture of 10"' K. The total luminosity from an HII region 
of radius Rq is then proportional to the volume of the HII 
region with a correction factor to take into account the 
other ionization stages of carbon. 

Because we deal with a fine, the opacity may not be 
small depending on the gas distribution and viewing ge- 
ometry. In fact, opacity effects are important at large in- 
clination angles. To obtain an edge-on view of the model 
galaxy in the ^-F3/2 — ^ P1/2 C"*" line, we have made an ac- 
curate calculation of the radiative transfer in this line. For 
each line of sight through the disk, we sample the line pro- 
file with bins of 1 kms~^ width, and calculate the emergent 



intensity in each velocity bin, including saturation effects. 
We assume that the intrinsic velocitjf dispersion of a PDR 
is 1 kms^^. For the diffuse medium, the velocity dispersion 
is 10 kms^^. We do not account for absorption in the far 
infrared continuum. We obtain tc+ = 0.40 for an edge-on 
view, and Tc+ = 0.10 for a face-on view. 

3. 7. Implementation of the model 

To avoid transient stages of the simulation, the code is 
evolved during a few time steps. We stop the simulation 
when stable results are obtained on a time scale of 20 
Myrs. This time scale corresponds to about half the life- 
time of a giant molecular cloud before disruption by pho- 
toevaporation. This is the reason why we can not integrate 
further in time without treating gas recycling. The time 
step has been fixed at 10^ years, shorter than the lifetime 
of the most massive stars. We have checked the reliability 
of the calculations by different tests: 

- We have verified that the total UV luminosity from the 
stellar population emerges from the galaxy either at the 
same wavelength, or at far infrared wavelengths for the 
light reprocessed by dust grains. The total luminosity is 
conserved with an accuracy of 1 %. 

- When the number of pseudo-photons leaving each OB 
association n-y is too small, the map of UV radiation field 
is noisy with a few extremely bright spots and large voids. 
This is due to undersamphng of the galaxy volume. The 
number of pseudo-photons should be as large as possible, 
but we have verified that we obtain a good map of the 
UV radiation field with 100 pseudo-photons per OB as- 
sociation. The map is smooth in the vicinity of the OB 
associations, hence the ratio of FIR emissions from the 
diffuse and dense gas stays constant with increasing n..^. 

- The cell size is also a critical parameter: since xuv 
is proportional to r^J^ it might be underestimated for 
smah clouds very close to OB associations. This has severe 
consequences for the C+ emission, which scales roughly 

''"leU^'^diXuv), but little or none for the FIR emission 
which varies as r'l^nXuv since in that case there is no 
resultant scaling with rceii- To test the validity of the 
adopted resolution, we have performed a run restricted to 
one quadrant only, with a cell size of 6.1 pc. We observed 
no large variation in the C"*" emission and thus conclude 
that the adopted resolution of 12 pc is correct for our pur- 
pose. Note that the volume filling factor decreases to 1.3 
% in the high resolution run, because we fill the space in 
a more accurate way using a higher spatial resolution. 

4. Results 

Table 2 summarizes the input parameters for the Stan- 
dard Model, and Table 3 presents the results. The star 
formation rate from 2 to 60 is fixed at 4.0 MQyr~^, 
with a star formation efficiency e of 5%, as defined in Sect. 
3.4. With these values, the modehed Ha luminosity and 
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Ha radial profile are in quite good agreement with the 
observed data (Kennicutt 1989). This is also true for the 
UV luminosity at 2000 A and the 6 cm luminosity. We 
are thus confident that the massive star population is well 
constrained by the observed data. With a low mass cut- 
ofi^ at 10 Mq stars, we overestimate the UV luminosity 
at 2000 A of the modelled stellar population by 20%, be- 
cause wc miss the contribution to the UV continuum of 
lower mass stars, between 2 and 10 Mq. 

4-1. Disk opacity in the UV 

We have computed an average opacity over the galaxy in 
the UV and for the Ha line. We define this opacity as : 

r = -ln(U'"'erge„t liemitted-^ ^j^g^g ^emitted ^j^g ^Q^g^j 

luminosity in the disk at a given wavelength and £<^™<^r9ent 
is the emergent luminosity. This opacity is computed for 
two different viewing angles of the model, i = 0° for face- 
on and i = 90° for edge-on. We have found a significant 
opacity for the face-on view, at 1000 A, 2000 A and for Ha, 
namely r(lOOOA) = 0.8, t(2000A) = 0.7, T{Ha) = 0.60 
for the whole galaxy. 

The opacity is controlled simultaneously by the geom- 
etry of the molecular cloud ensemble and by the diffuse 
medium. If we ignore the extinction due to the diffuse com- 
ponent, we find an opacity of 0.51 at 1000 A. This value 
is due to geometrical effects, mostly blocking of the UV 
radiation by molecular clouds, and it does not depend on 
wavelength. Thus we can write the opacity at any wave- 
length in the UV as t\ = 0.51 -I- , the second term 
accounting for the wavelength dependence of the extinc- 
tion in the diffuse medium. 

A global opacity of r ~ 0.8 corresponds to a fraction of 
approximatively 45% of the far UV stellar radiation leav- 
ing the galaxy disk, mostly above or below the main plane. 
Most of these photons have not been scattered because the 
probability of leaving the disc after a scattering event is 
low. This significant fraction of the radiation from mas- 
sive stars leaking out of HE regions could contribute to 
the maintenance of the Reynolds layer of ionized gas. The 
derived face-on opacity at 2000 A, 0.7, falls well within 
the range of opacities derived by Buat & Xu (1996). The 
mean extinction in their sample of nearby spiral galaxies 
is ~ 0.9 mag at 2000 A. Though the opacity is not very 
large, the mean distance travelled by a UV photon before 
absorption is quite small, 440 pc, roughly equal to the HI 
disc thickness. As shown on Fig. 6, there are however pho- 
tons travelhng to much larger distances, 1 to 2 kpc, with 
a small probability (0.01). Conversely, many zones in the 
interarm receive very few UV photons. Due to the lower 
gas density, few OB associations are created in the inter- 
arm region. The numerous OB associations in the arms 
are too distant to contribute to the local radiation field 
since the arm/interarm separation is larger than 1 kpc in 
the disk. 



The distribution of xuv values provide further infor- 
mation on the radiation field resulting from the OB asso- 
ciations (Fig. 7) . Whereas most of the galaxy is exposed 
to a low UV radiation field, it is possible to find regions 
with high UV intensity {xuv > 1000) even at a moderate 
spatial resolution. The total dynamical range of the UV 
radiation field extends over more than 4 orders of magni- 
tude. This huge variation can be explained hy the close 
association of OB associations and molecular clouds: in a 
galaxy with a prominent spiral structure, OB associations 
are born in the spiral arms, where the gas density is the 
highest. This maximizes both the illumination of molec- 
ular clouds by UV radiation and the absorption of UV 
radiation by molecular gas, hence the heating of molecu- 
lar gas. For the model galaxy, we find that 30 % of the 
total number of cells with molecular gas are exposed to a 
strong or median radiation field (xuv > 10)- These cells 
are located in 40% of the molecular clouds. This figure is 
comparable to the clouds in Milky Way: Solomon et al. 
(1985) found that in the Galaxjf, at a resolution greater 
than 10 pc, 25 % of the molecular clouds are warm and 
associated with HII regions. Also, Williams and Mc Kee 
(1997) estimate that at least one OB star is present in half 
of the giant molecular clouds with masses larger than 10^ 
Mq . The probability to find massive stars or clusters asso- 
ciated with a giant molecular cloud increases sharply with 
the cloud mass and reaches almost 1 for masses larger than 
8 X lO'^ Mg (Wilhams & MacKee 1997). Our numerical 
results are in agreement with these facts. 

4-2. Far infrared emission 

We now discuss the emergent radiation from the model 
galaxy and start with the infrared emission. As in the 
Desert et al. (1990) dust model, the luminosities in the 
IRAS bands are computed as AnD'^vSi,, where Sjy is the 
total observed fiux density and D is the distance to the 
object. The infared colors are given as the fiux density 
ratios, to compare with observed data. 

The model galaxy has very similar emissions as 
NGC 6946 at 60-100 & 200 //m, with outputs of 5.1, 8.5 
and 4.9 x lO'' Lg, corresponding to 114%, 128% and 144% 
of the luminosities observed at those wavelengths. The 
far infrared emission comes from both the molecular and 
atomic gas phases. 

The UV radiation is the main heating mechanism of 
the dense and diffuse gas phases, with contributions of 4.2 
X lO*' Lq and 6.1 x 10^ Lq at 60 & 100 ^im. The contri- 
bution to the FIR emission of the old stellar population, 
described by the ISRF, is a factor 3 lower, with 0.9 x 10^ 
Lq and 2.4 x 10^ Lq in the 60 & 100 bands. The sit- 
uation is different at 200 /xm, where dust grains heated by 
the UV radiation or by the ISRF have almost equal contri- 
butions to the total luminosity: 3.0 x 10^ Lq for the UV 
and 1.9 X 10® Lq for the ISRF. The contribution from the 
inner parts of clouds illuminated by the attenuated ISRF 
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is only 0.5 x 10^ Lq. As a whole, 72 % of the far infrared 
luminosity can be attributed to UV heated gas, which is 
mostly molecular. The remaining 28% corresponds to dust 
heated by the ISRF, at locations far away from the OB 
associations. 

The diffuse and dense phases have similar contribu- 
tions to the total FIR emission, with a shght excess from 
the molecular clouds, 54 % versus 46 % from the diffuse 
gas. This significant contribution from the diffuse gas is 
due to the fact that it occupies a large fraction of the 
galaxy volume. Hippclcin ct al. (1996) also conclude from 
ISO observations of other nearby galaxies (M51, MlOl) 
that the neutral atomic gas has an important contribu- 
tion to the far infrared emission. The contribution from 
the atomic gas may be underestimated because we do not 
take into account the compression of the diffuse gas in the 
spiral structure. Comparing with molecular clouds, we can 
estimate that, having atomic gas concentrated in the spi- 
ral arms would result in a brighter FIR emission, with 
a slighty warmer color temperature since the dust grains 
would be closer (in average) to the heating sources. A pre- 
cise estimate of the magnitude of the effect is beyond the 
scope of this paper. 

The global infrared excess for the model galaxy, IRE, 
is defined as the luminosity ratio IRE= Li2-iooiJ,m/ LLy, 
with L^y = N^ychvLy and /f./y£j,=13.6 cV. At the disk 
scale, the IRE takes the value 5.9, in agreement with ob- 
servations of Galactic HII regions (Caux et al. 1985, Myers 
et al. 1986). 

The diffuse and dense gas (atomic and molecular) have 
the following contributions to the total luminosity of the 
C"*" 158 /zm fine: 77% from the dense phase and 23% from 
the diffuse phase. Less than 10** Lq comes from HII re- 
gions. The total emission of the galaxy is 2.5 x 10^ L©, 
a factor 2.5 lower than the measured value, 6.3 x lO'^ 
L© (Madden et al. 1993). Compared to the 60-100 /im 
far infrared emission, the 0+ fine represents 0.21% of the 
FIR (60-100 //m) emission. This figure is comparable to 
the observed ratio for other spiral galaxies with 0.1 - 1 % 
(Lord et al. 1996). Nevertheless, the value for NGC 6946 
was found to be 0.6 % (Madden et al. 1993), and in the 
Galaxy, Shibai et al. (1991) and Wright et al. (1991) have 
measured Lc+ = 0.7 % Lpjji with the same definition of 
^FiR as above. 



4-3. Radial profiles 

The 60 /im radial profile is shown on Fig. 8a. There is a 
large decrease from the inner to the outer parts of the disk, 
about two orders of magnitude. In NGC 6946, the same 
behaviour has been observed by Tuffs et al. (1996) using 
ISO. Averaged over the model, the Seo/Sioo infrared color 
appears to be slightly different in the two phases: 0.32 
for the diffuse phase and 0.38 for the dense phase. This 
FIR color decreases with increasing radius from 0.40 in 



the center to 0.23 at R ~5 kpc (Fig 8b), in agreement 
with the maps by Engargiola (1991). The decrease is seen 
in both phases, with Seo/Sioo ranging from 0.35 to 0.23 
for the diffuse phase, and from 0.42 to 0.28 for the dense 
phase. 

The radial profile of the intensity of the 0+ ^^3/2 - 
158/xm line (Fig. 8c) shows a much flatter gradient 
than the FIR emission. This is due to the logarithmic de- 
pendence of the line intensity on the incident radiation 
field in PDRs. As shown in Fig. 8c and 8d, the diffuse 
atomic gas is the main source of C+ emission at large dis- 
tance from the nucleus, for radii larger than 4 kpc. It is 
thus possible to determine the intrinsic Lc+ /I^fir lumi- 
nosity ratio from the two gas phases, using the data at 
R ~ 2kpc for the molecular gas and data at R> 5kpc for 
the atomic gas. We find that Lc+fl^FiR is equal to 0.10 
% in the diffuse phase and to 0.25% in the dense phase. 



44- Maps 

We show on Fig. 9 face-on maps of 100 /im, C+ and 
UV(912-2000 A) emissions. Edge-on maps at the same 
wavelengths are shown in Fig. 10 for comparison. Com- 
pared to the 0+ observations of the edge-on galaxy NGC 
891 (Madden et al. 1994), there is an overall agreement. 
In particular, the scale height in C+ is predicted to be 
larger than the scale height of the CO emission, due to 
the contribution of the diffuse neutral and ionized media 
which have a larger scale height (Fig. 11). 

In the face-on 0+ map, there is a large hole in the inter- 
arm regions in the NW, at a similar position to the hole de- 
tected by Madden et al. (1993) with the KAO. This hole is 
due to the lower density of molecular gas and of OB associ- 
ations in the interarm regions. Therefore few UV photons 
illuminate this region and the radiation field is very low. 
The map shows many details and a large contrast between 
arm and interarm regions. We have smoothed the image 
from the model to the resolution of the KAO observations 
(50" beam = 1.2 kpc at the distance of NGC 6946). The 
contrast between the brightest regions and the disk drops 
by a large factor. This resolution effect may explain the 
low dynamical range found in the observed data. If PDRs 
are the main source of C"*" 158/iTO radiation in galaxies, 
we predict that the emission should have more contrast 
at higher spatial resolution. This could be tested by maps 
of external galaxies made with the future Stratospheric 
Observatory For Infrared Astronomy (SOFIA). 

The edge-on maps at 100 /im and in C+ are fairly 
symmetrical with respect to the center. The edge-on 0+ 
map shows however a hole in the central region ( r < 
500pc) which does not appear on the 100 jim map. This 
hole is largely due to the large opacity for these lines of 
sight (tc+ = 0.4). 
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4-5. Sensitivity of the model to input parameters 

The model results are of course sensitive to the input pa- 
rameters, therefore we have run different models deviating 
from the standard model by one parameter. 

Because of the poor knowledge of the albedo in UV, we 
have run a model with a lower albedo of dust grains, uj= 
0.4. We find that the opacity increases to 1.0 at 1000 A & 
0.90 at 2000 A. The 60 /im emission from the dense phase 
increases by 5%, while the 100 and 200 fim emissions both 
decrease by 10%. This difference in far infrared emission 
is due to the moderate increase of the opacity which leads 
to a warmer dust temperature. The effect on the emission 
from the diffuse phase is negligible. 

A more extreme case is for a null albedo, suppressing 
any scattering effect. In that case, wc maximize the UV 
opacity and the FIR emission. The opacity increases to 
1.25 at 1000 A and 1.01 at 2000 A respectively. As a result 
of this larger absorption, the 60-200/um emission increases 
by 47 % as compared to the standard model. 

In another run, we have kept the total mass of molec- 
ular gas constant, but used a lower mean density, 20 H2 
cm~^ instead of 50 H2 cm~^, to increase the clouds sizes. 
The volume filling factor is then 3.8 %. These larger clouds 
block more Hght, and 30 % only of the molecular cells are 
heated, instead of 40% in the standard model. As a result, 
the 60-200 (im luminosity decreases by 10%, to 16.7 x 10^ 

We have also investigated the effect of the number of 
OB associations: we have kept the same star formation 
rate but have gathered adjacent associations to form more 
powerful sources. As a consequence, noB decreases from 
12000 to 3000. Then a smaller fraction of the cloud popu- 
lation is heated, 15%, as compared to 40% in the standard 
model. But because these cells are heated by more power- 
ful OB associations, the far-infrared emission is larger and 
reaches 20.4 x 10® Lq. Thus the FIR emission depends 
slightly on the number of associations. The C"*" emission 
decreases to 1.8 x 10'' L©, because of the smaller number 
of illuminated clouds. 

If we now increase the SF efficiency, from 5 to 10%, 
so as to double the UV luminosity, the production rate 
of Lyman continuum photons increases by 80 %. In that 
case, the mean UV opacity is 0.78. The FIR 60-200 nm 
luminosity increases by 55 % to 28.7 x 10® Lq. This shows 
that the FIR emission is not a linear function of the UV 
luminosity in our model. This non-linear behaviour arises 
because the opacity is largely controlled by geometrical 
effects. With a larger star formation activity, HII regions 
are very large and can destroy molecular clouds efficiently. 
Thus the mass of molecular gas decreases in the model 
with a higher SFR. This is the main reason for the non- 
linear behaviour. This result has been established with the 
same number of OB associations, while an increased SFR 
win probably lead to more associations in the disk. How- 



ever we have previously shown that the FIR emission does 
not depend strongly on the number of OB associations. 

We have investigated the effect of the atomic density 
on the size of HII regions, because we probably overes- 
timate the diameter of HII regions, using a mean atomic 
density and neglecting the dust absorption. If the local gas 
density is multiplied by two, the volume of the Stromgren 
sphere is 4 times smaller than in the standard model. The 
60-200 /im luminosity of dense molecular gas increases by 
10% to 1.1 X 10^° Lq. This is explained by the reduced 
destructive effect of HII regions on molecular clouds, and 
then the larger chance for photons to be absorbed by 
molecular gas. The respective contributions from the dif- 
fuse and dense gas to the FIR(60-200 fim) are now 34% 
and 66%. 

This last test shows that the distance between clouds 
and OB associations has a strong infiuence on the UV re- 
processing by dense gas. For the standard model, we have 
calculated the mean distance between an OB association 
and the nearest cloud edge, doB /clouds and have found a 
value of 35 pc, the mean distance between clouds centers 
is 37 pc. To have a larger separation between clouds and 
OB associations, we have increased Vescape to 30 kms~^. 
We obtain doB/cioud= 39 pc. The FIR(60-200 fxm) emis- 
sion from the dense phase decreases by 15% because of 
the smaller solid angles of the clouds viewed from the as- 
sociations. As for the FIR (60-200 fim) from the diffuse 
phase, it slightly increases by 6%. 

We have shown that part of the UV opacity is due to 
geometrical effects. Indeed the UV opacity is lower when 
the molecular clouds are distributed uniformly in the disk. 
We have used an earfier epoch of the simulation, when the 
distribution of gas clouds is axisymmetric. We have kept 
the same value for the other parameters (number of OB 
associations, star formation rate, etc.). In that case, the 
clouds occupy a larger fraction volume of the galactic disk, 
and the mean distance between clouds increases. Because 
of this larger mean distance between clouds, the opacity 
at 1000 A decreases to 0.47. 

5. Discussion &; conclusions 

We have shown that with simple assumptions about the 
birth of massive stars and their relationships with the 
ISM, we can reproduce qualitatively and quantitatively 
the characteristics of the UV, Ha and FIR emissions of 
a particular object.the Sc galaxy NGC 6946. For such a 
galaxy with a prominent spiral structure, having a large 
mass of neutral gas, and forming stars actively, the ob- 
served far infrared emission is produced both in molecular 
gas and in the diffuse atomic gas. More precisely, 54 % of 
the FIR (60-200 //m) emission comes from dust grains in 
giant molecular clouds. Dust in the diffuse neutral atomic 
gas contributes to about 46 % of the total FIR luminosity. 

We have evaluated the respective contributions of the 
UV radiation from massive stars and of the radiation field 
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from the old stellar population. We find that 72 % of 
the FIR luminosity can be attributed to UV heated dust 
grains, which reside mostly in molecular clouds envelopes. 
The remaining 28 % is due to dust heated by the radiation 
field from the old stellar population at locations far away 
from OB associations. 

We have calculated the emission of the model galaxy 
in the ^P3/2 ~^ -Pi/2 fine structure line of 0+ at 158 //m. 
In the spiral arms, photon dissociation regions at the sur- 
faces of molecular clouds are the main source of the emis- 
sion. We have found a large arm-interarm contrast in this 
line. This effect could be tested by high angular resolution 
maps of galaxies. It results naturally from the combina- 
tion of a lower gas density and lower radiation field in the 
interarm regions, because of the short mean free path for 
UV photons, ~ 440 pc. As a whole, PDRs represent 76% 
of the emission. The contribution from the diffuse phase 
is found to be ~ 24 %. Our model is able to account for 
about 40 % of the observed C+ emission of NGC 6946. 
The emission from PDRs should be viewed as a lower limit 
since we use the model by Le Bourlot et al. (1993) with 
low abundances of carbon and other elements in the gas 
phase: [C]/[H]= 3 x 10"^ The average value is 1.3x 10"" 
for Galactic diffuse clouds (Snow & Witt 1996), a factor 
of 4 larger than the value used in the model. Since the C+ 
158 iim emission scales roughly with the column density, 
hence the carbon abundance, the total C+ luminosity from 
PDRs could be larger by at least a factor three than our 
current model prediction. This would increase the contri- 
bution from PDRs to the total C+ emission of the model 
galaxy: with this scaHng factor, the predicted 0+ lumi- 
nosity of PDRs would reach 6 x \^ Lq. Moreover, the 0+ 
emission from the diffuse gas is overestimated, because the 
model of Wolfire et al. (1995) assumes [C]/[H] ~ 3 x 10"" 
in the gas. Thus the diffuse emission could be 2-3 times 
smaller than in our standard model. Adopting [C]/[H]=1.3 
X 10"'' in both phases would thus enhance the differences 
of Lc+ 1 LpiR between dense and diffuse gas. 

The knowledge of the cloudy nature of the ISM, and 
of the global structure of the galaxy, is important to de- 
termine how far UV photons can travel away from OB 
associations. The filling factor and the mass/radius scal- 
ing law appear to be major parameters for the transfer 
of stellar radiation in the galaxy disk, because thejf de- 
termine at the same time the obscuration and the size of 
the emitting regions. Other important parameters are the 
number of OB associations and the sizes of HII regions, 
because with a large number of OB associations or with 
small HII regions, molecular clouds are on average closer 
to massive stars, and are thus more efficiently heated. 

In all the models we ran, we have found that the aver- 
age internal UV opacity is of the order 0.8. The discs are 
therefore moderately opaque in the UV, as measured by 
Buat & Xu (1996). This moderate opacity holds for face- 
on discs. Edge-on discs are quite opaque, with a small 
fraction of the luminosity escaping, less than 1.0% of the 



face-on luminosity. This fraction corresponds to an equiv- 
alent extinction of 5 magnitudes in the UV. 

These results have been obtained using a crude de- 
scription of the interstellar medium. The adopted spatial 
resolution results from a compromise between astrophysi- 
cal requirements and computational needs, but is certainly 
very poor compared to the complexity of the interstellar 
medium. The good agreement of the observed and pre- 
dicted large scale properties shows nevertheless that the 
transfer of UV radiation, and the role of the radiation for 
the gas and dust heating, are correctly described at the 12 
pc scale. This is in agreement with previous works estimat- 
ing that dust heating by UV radiation occurs principally 
at large distances from massive stars (Murthy et al. 1992, 
Leisawitz & Hauser 1988). 
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Figure captions 

Fig. 1. C0(J=2-1) map at 13" resolution obtained with 
the IRAM 30m radiotelescope. 

Fig. 2. V-I image, obtained with the 1.20m telescope at 
the Observatoire de Haute Provence (P. Boisse, private 
communication). The field of view is 10' by 10', the pixel 
size is 2.3". The gray scale runs from white for blue colors 
to black is for red colors. 

Fig. 3a. R band image of NGC 6946 used as input for the 
calculation of the gravitational potential. Field stars have 
been removed, the image has then been rotated and de- 
projected to get a face-on view of NGC 6946. 
Fig. 3b. A face-on view of the model galaxy, with molec- 
ular clouds drawn as circles and OB associations drawn 
as stars. Only 20 % of the OB associations and molecular 
clouds are drawn. 

Fig. 3c. A close-up view of Fig. 3b. The clouds are drawn 
at their exact size in the model. 

Fig. 4. Adopted rotation curve (dot-dashed line), angular 
frequency fi(r) (full line), and 0, ± ^ curves. The corota- 
tion for the adopted pattern speed is located at r = 3.5 
kpc. 

Fig. 5. Distribution of the production rate of Lyman con- 
tinuum photons N^yc, for the modelled population of OB 
associations (full line) in NGC 6946, and for the Galaxy 
(dot-dashed curve) (McKee & WiUiams (1997)). The plot 
can be read for instance as 500 OB associations out of 
12000 have ^Lyc > s-\ 

Fig. 6. Distribution of the distance from their parent OB 
association, travelled by UV photons before absorption. 
The median distance is 120 pc, and the mean distance is 
440 pc. 

Fig. 7. Distribution of the UV intensity measured rela- 
tive to the ISRF in the modelled galaxy, for cells filled 
with diffuse gas (dot-dashed line) and for cells filled with 
molecular gas (full line). 

Fig. 8a. Radial profiles of the 60 fim surface brightness. 
The thin solid line shows the combination of the two gas 
phases, while the triangles show the contribution of the 
diffuse atomic gas, and the squares the contribution from 
the molecular gas. The thick solid line shows the mean 
radial profile deduced by Engargiola (1991) from an IRAS 
image. 

Fig. 8b. Radial profiles for the ratio of 60 and 100 /xm 
fiuxes, global and for the two phases. 
Fig. 8c. Radial profiles of the C"*" emission: global and for 
the two phases. 

Fig. 8d. Radial profiles of Lc/L(6o-ioo^tm) ^ global and for 

the two phases. 

Fig. 9. Face-on views of: 

a) 100 fxm emission, at 48 pc resolution. The gray scale 
ranges from 1 to 10^ Lpj pc^^. We have overlaid contours 
of the same map convolved with a 750 pc beam: the levels 
are at 10, 30, 60, 100, 200 Lq pc'^ 

b) C+ fine, at 48 pc resolution. The gray scale ranges from 



10~^ to 1 L© pc~^. Overlaid contour levels from 0.1 to 0.6 
by 0.1 L© pc^^ for the same image convolved with a 750 
pc beam. 

c) emergent UV surface brightness (resolution: 48 pc). The 
gray scale ranges from 10~^ to 10^ L© pc~^. 
Fig. 10. Edge-on views of NGC 6946, at 48 pc resolution. 
The linear scale is not identical for both axes. 

a) 100 jjm IRAS band, with contour levels at 10^, 3 x 10^, 
6 X 10^ and lO'' L© pc^^. 

b) 0+ , at the same resolution, accounting for the fine opac- 
ity. The gray scale ranges from to 50 L© pc~^. Contour 
levels from 10 to 50 by 10 L© pc~^. 

Fig. 11. Average vertical profiles through the disk for the 
C"*" (dot-dashed line) and C0(l-O)(full line) emissions. 
The CO (1-0) profile presents a smaller scaleheight (40 pc) 
than the C"*" profile (100 pc). We have assumed that the 
CO(l-O) emission is proportional to the molecular gas col- 
umn density. 
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Table 1. Observed parameters for NGC 6946. 



Observed data 


Value 


Reference 


Distance 


5 Mpc 


De Vaucouleurs (1979) 


Inclination 


34° 


Considere & Athanassoula (1988) 


Mh2 (R < 6') 


2.0 10" Mq ^ 


Young & Scoville (1982) 


Mhi (R < 6') 


2.6 10'' Mq ^ 


Boulanger & Viallefond (1992) 


Ll950-2050A 


5.4 10** Lq ^ 


Buat et al. (1989) 




1.0 10* Lq 


DeGioa et al. (1984)) 


B? 


8.49 mag 


RC3 


IRAS Ufim 


2.3 lO'' Lq ± 20 % ^ 


Engargiola (1991) 


IRAS 25/im 


1.4 10'' Lq ± 20 % = 


Engargiola (1991) 


IRAS 60/um 


4.5 lO'' Lq ± 20 % ' 


Engargiola (1991) 


IRAS 100/im 


6.6 10" Lq ± 20 % ^ 


Engargiola (1991) 


IRAS 200/im 


3.4 10" Lq ± 20 % ^ 


Engargiola (1991) 


Lc+ 


6.3 10^ Lq " 


Madden et al. (1993) 


Qthermal 
^6cm 


83 mJy ± 25 


Klein et al. (1982) 



We use the value 3.8 x 10^^ erg s ^ for the solar luminosity at any wavelength. 



^ We use Nh2/Ico(i-o) = 2.3 10^° mol cm~^/(K kms~^) and a total intensity Ico(i-o)= 569 Kkms~^ in a 45" beam, we do 
not account for projection effects. 

^ We use Mhi = 1.9 x 10^° Mq at a distance of 10 Mpc, corresponding to 4.7 x lO" Mq for a distance of 5 Mpc. By performing 
an integration over the HI map, we estimate that the part with R < 6' contributes to 55 % of the global emission. 
^ The 2000 A flux, 6.92 x 10"^^ erg cm ^s ^, is corrected for the Galactic extinction and integrated over a 100 A band. 
Note that Donas & Deharveng (1987) have reported a flux at 2000 A twenty times lower (1.66 x 10"^^ erg cm^^s^^A"^). 
Measurements at other UV wavelengths are available in Rifatto et al. (1995). 

* We use the uncorrected flux iHa= 1.31 x 10^^" erg cm^^s^^ and the relationship Lhc = 3.13 x lO^'' ^Mpc fea from Young 
et al. (1996). We have subtracted the contribution from the nucleus, estimated to be 20 % of the total luminosity. Note that 
Young et al. (1996) give a lower total flux, 3.38 x 10~^^ erg cm~^s~^, while Kennicutt (1989) reports a mean surface brightness 
of 4.2 X 10^^ erg s^^ pc^^ corresponding to a total Ha luminosity of 3 x lO'^ Lq at the adopted distance. 

® We estimate the luminosity in the band centered on frequency by i^F^. We use the IRAS and KAO fluxes reported by 
Engargiola (1991) for a radius R < 5.6' and correct from the contribution of the nucleus (central 45") to obtain 12 Jy at 12/im, 
15 Jy at 25/im, 114 Jy at 60/im, 283 Jy at 100/im and 288 Jy at 200/im. The global color for the disk is Seo/Sioo = 0.40. The 
far infrared luminosity is Lfir( 60-100^^)= Leo^ + Lioo/x = 1.1 x 10^° Lq, hpjjK^ 60-200^^) = 1-45 x 10^° L© 
^ We have excluded the contribution of the nucleus and rescaled the data by Madden et al. (1993) for the adopted distance. 

We use an estimation of 625 mJy for the total flux at 6cm (Klein et al. 1982), subtract the nuclear contribution estimated to 
be 33% of the total flux, and we keep 20% of the resultant disk flux as the thermal component. 



Table 2. Input parameters for the standard model. 



Parameter 


Adopted value 


cell resolution 


12.2 pc 


HI half thickness 


180 pc 


H2 scaleheight 


65 pc 


mean Vos 


10 kms"^ 


n-y 


100 


time step 


10^ years 


number of molecular clouds 


18000 


number of OB associations 


12000 


star formation rate (2 - 60 Mq) 


4.0 Mq yr"^ 


volume filling factor of H2 


2.5 % 
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Table 3. Output values from the simulation for the standard model. 



Parameter 


Total 


Dense 


Diffuse 


uv 


ISRF 


L912-2000A (L©) 


1.3 10'" 


— 


— 


— 


— 


Ll930-2030A (Lq) 


5.2 10** 


— 


— 


— 


— 


L//,> (Lo) 


1.2 10® 


— 


— 


— 


— 


B° (mag) 


8.96 


— 


— 


— 


— 


L12 (Lq) 


2.1 10" 


1.0 10" 


1.0 10" 


1.7 10" 


0.3 10" 


L25 (Lq) 


2.0 10" 


1.0 10" 


1.0 10" 


1.7 10" 


0.3 10" 


Leo (Lq) 


5.1 10" 


2.9 10" 


2.2 10" 


4.2 10" 


0.9 10" 


Lloo (Lq) 


8.5 10" 


4.5 10" 


4.0 10" 


6.1 10" 


2.4 10" 


L200 (Lq) 


4.9 10" 


2.5 10" 


2.5 10" 


3.0 10" 


1.9 10" 


L60-200 (Lq) 


18.6 10" 


9.9 10" 


8.6 10" 


13.3 10" 


5.2 10" 


Lc+ (Lq) 


2.6 lO'' 


2.0 lO'' 


0.6 lO'' 






Seem (mjy) 


124 










Lc+/Lf/r (eo-ioo) 


0.18 % 


0.25% 


0.10 % 







All luminosities are emergent luminosities. The total UV luminosity generated by 
OB associations in the disk between 912 and 2000 A is 3.10 x lO" Lq. 
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